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Abstract 


The paper presents the possibility of stabilizing higher harmonics of piston speed generated by torsional 
vibrations of the shaft. As a result of torsional vibrations, the transient speed of the piston deviates from the values 
resulting from the known formulae describing geometrical dependences of the piston position versus the angle of 
rotation of the crankshaft [1]. This phenomenon causes the desired effect of damping torsional vibrations. As it is 
known, in larger engines of the cylinder diameter exceeding 120mm and the number of in-line cylinders greater than 
or equal to 6, structural damping is insufficient and it becomes necessary to use torsional vibration dampers. In 
article [2] attention was drawn to the effect of damping higher harmonics of the moment generated by the engine 
whose lateral surfaces are coated with layers of nanotubes, which was not however, the main subject of the article. 
This work presents a preliminary analysis and determinants of the efficiency of vibration damping due to the impact of 
nanotubes on the reduction of the amplitude of higher harmonics of the moment generated by the engine. In addition, 
the possible mechanism of the phenomenon of vibration damping by a layer of nanotubes is presented. However, one 
should emphasize that at this stage the authors do not close the debate concerning the mechanical properties of the 
structures based on carbon nanotubes (CNTs), but in fact, they open up such a discussion. Besides the confirmed 
properties of nanotubes, such as extremely high tensile strength of properly structured fiber formed by CNTs, there is 
very little data concerning the properties of chaotic structures in which nanotubes combine in larger structures 
adopting any directions. No less important is the base on which nanotubes are grown, and it can be very different. The 
authors of the paper have presented preliminary results of a positive impact of selected properties of CNTs on the 
vibrations of the crankshaft. 


Keywords: friction in internal combustion engine, nanotubes layers, crankshaft torsional vibration 


1. Introduction 


One of the methods of reducing friction losses in a combustion engine is coating the lateral 
surfaces of pistons with proper layers with a low coefficient of friction. This method assumes that 
in a combustion engine the lateral surface of the piston comes into direct contact with the cylinder 
bearing surface. As it is confirmed by the research and computer simulations, the described 
situation is a phenomenon of a relatively short duration, usually when the engine is started, 
following a longer standstill [3]. Consequently, in a properly designed piston-cylinder group, the 
friction losses depend to a small extent on coating the lateral surfaces of the piston with layers 
with a low coefficient of friction. The importance of these coatings 1s exposed in easier engine 
start, and first of all, in reducing abrasive wear, before a continuous oil film is formed on the 
cylinder bearing surface. However, it appears that coating the cylinder bearing surface with 
specific layers may lead to reducing vibrations of the crankshaft, and thus the whole sequence of 
drive transmission. Fig. 1 presents schematically and disproportionately the combination of the 
side-surface of a piston coated with a layer of carbon nanotubes with the cylinder bearing surface. 
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Fig. 1. The combination of the side-surface of a piston coated with a layer of carbon nanotubes 
with the cylinder bearing surface — description in the text 


The structure of the lateral surface of piston 1 was outlined while taking into account only the 
trace of the curvature of the cutting tool that is left after finishing. Both on the lateral surface of the 
piston and on cylinder bearing surface 2 there is micro-roughness which is not relevant to the 
discussed issue, therefore, Fig. 1 pictures the plateau of these surfaces. In the grooves of the side- 
surface of the piston, there are embedded carbon nanotubes 3, which in fact form a chaotic 
structure. But in order to orientate the reader which direction the layer of nanotubes is deforming, 
they were pictured by horizontal lines. The layers of the oil film formed between the surfaces in 
contact was marked by vertical lines. The position of the piston on the left side in Fig. 1 


corresponds to the zero relative velocity of the piston and sleeve. The middle diagram shows 
laminar motion of the oil film layers when the piston is ascending, and the diagram on the right 
side describes the piston motion in the opposite direction. As a result of the laminar movement of 
the oil film layers, within the contact area of a layer of nanotubes, there occurs a tangential force 
deforming nanotubes located in the grooves of the lateral surface of the piston. The deformation of 
the nanotube layer may be elastic, elastic-plastic or plastic. The results of the test carried out by the 
authors of this article confirm that the structure of randomly oriented nanotubes mainly exhibits 
elastic properties. This means that the energy of deformation of the nanotube layer is absorbed at 
the very moment when the piston accelerates, and when the piston slows down its speed in relation 
to the sleeve, the energy is transferred onto the piston and further onto the power receiver. As a 
result, there occurs suppression of higher harmonics describing the movement of the piston, which 
is equivalent to damping torsional vibration of the shaft. As it is known, the basic component of 
damping torsional vibrations in a classic internal combustion engine is generated by the viscous 
force of the oil film on the cylinder bearing surface. Unfortunately, such damping of torsional 
vibrations causes large energy losses because the change in the velocity of oil layers which are 
moved laminarily causes energy absorption, whether the piston accelerates or slows down; 
therefore, irrespective of the direction of the piston motion as well. 


2. The course of engine's moment of resistance to motion 


The actual course of the coupling moment of an internal combustion piston engine with 
adriving machine deviates considerably from the calculation values based on geometrical 
relationships. The reason lies in numerous phenomena, but the vibrations of the following unit: the 
engine shaft — the driving machine shaft are the most important. From a very extensive study, we 
selected the most representative measurement results of the coupling moment of a two-cylinder 
combustion engine with a driving machine and the obtained results are presented in Fig. 2-4. 
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Fig. 2. The course of the coupling moment of a two-cylinder combustion engine with the machine driving the engine at 
the angular velocity of 1000 rpm for standard pistons whose lateral surface is not coated with a layer of nanotubes, 
description in the text 
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Fig. 3. The course of the coupling moment of a two-cylinder combustion engine with the machine driving the engine at 


the angular velocity of 1000 rpm after 10-minute grinding-in of nanotubes on the lateral surface of the piston 
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Fig. 4. The course of the coupling moment of a two-cylinder combustion engine with the machine driving the engine 


at the angular velocity of 1000 rpm after 60-minute grinding-in of nanotubes on the lateral surface of the piston 


The following are marked in the drawings: 

— the blue line — the course of the measured coupling moment, wherein the scale of the graph 
is described on the left vertical axis. The mean coupling moment is given on the right side 
of the drawing, next to the word “pomiar”, 

— the green line — the course of the processed signal of repeatable moment every two 
rotations of the crankshaft; the scale of the graph was described on the right vertical axis, 

— the black line — the course of the filtered signal causing unrepeatability of the coupling 
moment at the rhythm of processes taking place in the engine, the scale of the graph was 
also described on the right vertical axis. 

Due to the engine type and its modifications, the course of the coupling moment should in 


principle be repeatable every single rotation, but because of the difference in distance of two 


10 


double cranks from the driving machine there must occur certain differences in the generated 
moment of resistance related to the compression of the agent within cylinder volume. This results 
from differential susceptibility of drive transmission. 

What is important is the answer to the question of the cause of unrepeatability of the measured 
coupling moment. As one can notice, in the case of full repeatability of the coupling moment the 
filtered signal shown as the black line in Fig. 2-4 should amount to 0. However, since the filtered 
signal shows variable values, one should accept a definition of a certain indicator which would be 
the measure of unrepeatability of the course of the coupling moment. The simplest and, as it 
seems, the most proper indicator of the unrepeatability of the coupling moment is the sum of 
deviations from the mean value in relation to the number of included values. In this case, the mean 
value is 0, and the deviation is always assumed positive. The indicator of unrepeatability of the 
courses of the coupling moment is to be referred to as a in later parts of the article. 

As mentioned above, it is torsional vibrations that are the main cause of unrepeatability of 
the courses of the coupling moment. The frequency of these vibrations is dependent on the mass 
moments of inertia of the elements coupled with the shafts and the rigidity of the units connecting 
the shafts, and not on the angular speed of the set: combustion engine — driving machine. One can, 
therefore, put forward a hypothesis that the smaller the value of indicator a of unrepeatability of 
courses, the less the energy of torsional vibrations, and therefore less stress loading the engine's 
shaft. 

In Fig. 2-4 the values of indicators a were marked next to the word ,,Filtr” on the right side of 
the graphs. The value of indicator a is given in [Nm]. More details of the process of evaluating the 
value of the filtered signal are provided presented in literature [2]. 

In the case of a standard piston engine in which the piston lateral surface is not coated with a 
layer of nanotubes, the value of indicator a was 0.123 [Nm] — Fig. 2. Coating the lateral surface of 
the pistons with a 5 um layer of nanotubes brought about a decrease in the initial value of indicator 
a to the level of 0.113 Nm — Fig. 3. This value was obtained following a 10-minute running-in of 
the engine, since the moment when pistons coated with a layer of nanotubes were embedded. The 
shape and location of the layer of nanotubes is shown in Fig. 5. 


Fig. 5. The piston with a layer of nanotubes applied to the lateral surface prior to testing 


Further grinding-in the layer of nanotubes on the lateral surface of the piston during the period 
of 30 minutes caused a systematic decrease in the average value of the coupling moment, followed 
by its stabilization. The course of the coupling moment that was obtained after 60 minutes is 
shown in Fig. 4. The value of indicator a after 60 minutes of movement decreased to 0.097 [Nm] 
so in relation to the value a for standard pistons, the decrease in free vibrations of the set was 
reduced by approximately 25%. 
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Having completed the testing, a disassembly of the engine was performed. Examination of the 
pistons demonstrated that the layer of nanotubes was worn off to the edges of vertices that were 
left as a result of machining — Fig. 6. 


Fig 6. A section of the lateral surface of the piston following a 60min running process 


One should emphasize that applying a layer of nanotubes on the surface of an aluminum alloy is a 
process that is still in its experimental phase, and first attempts to control such a process gave 
positive results in very few laboratories [4-14]. Fig. 7-8 present the structures of the chosen ways 
of growing nanotubes on the surface of a carrier. 


Fig. 7. Pictures from a transmission electron microscope presenting CNTs 
with the embedded platinum nanoparticles on them 


Conducting the process of growth while applying a strong electric field allows for vertical 
growth of nanotubes, which is very important from the viewpoint of mechanical properties of the 
applied layers on the side-surfaces of pistons. The vertical growth of CNTs was presented in 
Fig.ures 8A), B) and C), whereas the anisotropic one in Fig. 8 D). 
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Fig. 8. Pictures from a scanning electron microscope (SEM) showing metal surface coated with CNTs 
with the embedded platinum nanoparticles on them 


Current research focuses on growing more durable layers of CNTs of the desired mechanical 
properties on the surface of aluminum alloys. 


3. Summary 


The research presented in the article shows the possibility of using coverings of layers of 
nanotubes on the lateral surface of the aluminum alloy piston in order to reduce unwanted torsional 
vibrations in the sequence of drive transmission. Further tests modifying the structures of 
nanotubes are being carried out in order to reduce frictional resistance and to more effectively 
reduce the vibrations of the engine shaft. The experiments concerning applying various elements 
and chemical compounds on CNTs suggest that damping properties of nanotube layers and the 
frictional forces generated by them on the contacting point of oil film allow for obtaining 
significantly better performance of piston combustion engines. Frictional losses of pistons on the 
cylinder bearing surface will, however, depend first of all on the parameters of the oil film that 15 
generated by properly designed lateral surface of the piston. Any layers applied to the lateral 
surfaces of pistons may reduce friction only in the phase of starting an engine, before oil film is 
formed on the cylinder bearing surface. 
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Abstract 


This study deals with the problem of the diagnostics of self-acting valves installed in reciprocating machines. 
Piston compressors and their valves are characterized. Compressor diagnostic systems, those offered and used, are 
described, including their restrictions. The signal generated by closing the valve is estimated. The authors have 
examined the relationship between the vibration signal and the condition of plate delivery valve of an air compressor. 
Simulations were made of the spring and plate wear and valve assembly errors. It has been shown that there is a 
relation between the selected vibration measure, mean peak-to-peak value of vibration accelerations and the technical 
condition of the discharge valve and measurement point. The authors have concluded, that in case of plate valves, one 
should expect that values of vibration measures depend on the measurement point and will decrease depending on its 
technical condition. A rule of diagnostic inference was formulated. 


Keywords: technical diagnostic, reciprocating machines, self-acting valves, condition of discharge valve, vibration 
signal 


1. Introduction 


Successful vibration-based diagnostics of rotor machines encourage manufacturers of 
diagnostic systems to develop and offer equipment and systems for piston (reciprocating) 
machines as well, including piston compressors. Various authors indicate possibilities of effective 
diagnostics of piston-crank mechanism components. However, certain problems occur when 
attempts are made to diagnose valves. 


2. Compressors and their valves 


There is a wide variety of design solutions in piston compressors. Also, various classifications 
exist. Depending on the piston-crank mechanism, there are compressors with and without crank 
mechanisms. Depending on the arrangement of cylinders (shape of the crankshaft) compressors 
can be of in-line, boxer, V-, star- and vane types. If we consider the shape of piston, there will be 
single-acting, double-acting, differential and stage-differential pistons. Due to piston guiding, 
compressors are with or without crossheads. Depending on the piston axis position, compressors 
may be horizontal or vertical. 

The medium and required pressure lead to this classification [3]: 

— non-critical use: low pressures, gases are not dangerous; 

— semi-critical use: high pressures, gases are not dangerous; 
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— critical use: media are potentially dangerous. 

Considering compressor dimensions, e.g. power N of compressor drive, compressors belong to 
these ranges: N « 50 kW, 50 kW< N « 100 kW, N 7 100 kW [3]. 

The medium being compressed results in another classification: air compressors, hydrogen 
compressors etc. 

For the control of compressor working process valves are used: self-acting and controlled 
(forced timing). Forced timing is used mainly for controlling the delivery rate of the compressor. 
The suction valve closing angle is controlled by an electronically controlled hydraulic cylinder. In 
most design solutions self-acting valves are used. 

Self-acting valves are of various types [4]: 

— plate valves; these valves, with a plate functioning as the sealing element made of metal or 
non-metallic materials; this type of valve is used in oil and gas, fertilizer and other 
industries; 

— ring valve; instead of a valve plate, independent rings are used. The rings are opened and 
closed by springs adjusted to the working conditions. The materials for rings are carbon 
fiber-reinforced composites, characterized by very high tightness, impact resistance, 
resistance to plastic deformation and particles, etc. Their significant advantage is that 
cylinder liner sliding surface is not damaged in case of ring fracture. Rings can be covered 
with an anti-adhesive material. Ring valves are intended for most extreme working 
conditions: hydrogen or light gases (molecular weight « 8 kg/kmol) in refineries and oil and 
gas industry; 

— mushroom valves; valves with a number of mushrooms in one assembly are used in 
compressors installed in natural gas transport systems where large quantities of gas are 
transferred at relatively low degree of compression, at low or medium flow rate, 180—600 
rpm. These parameters allow to apply high lift (max. lift: 8 mm); 

— concentric valve; used in single-acting compressors, concentric valves are the most effective 
solutions, comprising the whole piston surface. The suction and discharge valves are 
mounted in a joint body of cylinder head. Such solution is characterized by slight flow 
resistance and high durability. These valves are used in air and refrigerating compressors, 
sometimes in single-acting gas compressors; 

— lamellar valve. Lamellar (tongue) valves are generally an integral part of cylinder head. The 
lamella of the suction and discharge valve is integrated on one joint plate of the valve seat. 


3. Diagnostics of compressors 


Depending on the scope of identification of machine condition, condition monitoring and 
diagnostics are distinguished. As various equipment may be used, condition monitoring and 
diagnostics can be executed as remote condition monitoring and diagnostics, in-site monitoring 
and diagnostics or distributed condition monitoring or diagnostics. 

In practice [3], compressors are generally equipped with vibration transmitters sending signals 
to superior systems. TDC sensors are also installed in order to carry out measurements for various 
angles of the crankshaft by means of a portable measuring instrument. Vibration transmitters are 
usually mounted on the cylinder block. Major compressors are additionally equipped with 
vibration transmitters installed on the crosshead guide. Some compressors are fitted with indicator 
cocks, so that during their operation, if need arises, р-Г diagrams can be drawn by portable 
measuring instruments. During diagnostic inference the measured signals as a function of 
crankshaft rotation angle are analyzed. 

Those dealing with compressor diagnostics face a problem of identifying components of a 
vibration signal measured on the compressor cylinder. Normal vibration sources are impacts of the 
piston against cylinder liner and impacts of valves against the cylinder head, dependent on the 
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condition and loads. The difficulty lies in the fact that both impacts occur at a similar rotation 
angle of the crankshaft, Fig. 1. Other sources of vibration are possible. These activate themselves 
only after a defect occurs. In [3] a case is described where an increase of acceleration of cylinder 
vibration was considered as a symptom of incorrect valve condition, while inspection revealed that 
the cause was piston seizure in the cylinder. 

p Ta v 


T 


Fig. 1. P-v chart of a single-stage piston compressor: L; — work of the circulation, TDC — top dead centre, BDC — 
bottom dead centre (own study on the basis of [2]) 


The problem of piston impact on the cylinder liner has been discussed in [1]. The impact of 
self-acting valves requires a separate analysis. 


4. Examination of vibration signals from self-acting valves 


The impact of valve plate against the valve seat depends on the pressure difference before and 
after the valve, valve spring tension, plate lift and "elasticity" of the plate and the seat. The 
mentioned quantities may change due to valve wear. As a result, the impact force may: 

— increase, e.g. because of increased valve lift; 

— decrease, e.g. due to decreased spring stiffness. 

It should be expected that a diagnostic symptom of valves can be either an increase or decrease 
in the measure value of vibration generated by impacts in valves. The value of vibration measure 
also depends on the point of measurement. None of the authors on the subject indicates that the 
cylinder block is the best place for measuring vibration dependent on valve condition; points of 
measurement, however, are selected for their ease of mounting and conviction that the distance of 
vibration propagation of all possible sources is the shortest. 

A vibration signal from a plate delivery valve was examined, Fig. 2. 
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The measurements aimed at examining the influence of changes in the technical condition of a 
discharge valve on the second cylinder of a three-cylinder single-stage air compressor on the 
vibration signals of the valve, valve plate and cylinder head, Fig. 3. Vibration accelerations were 
measured. Vibrations on the cylinder block were not measured as its construction made it 
impossible to measure vibrations by means of available transmitters. 

Vibration accelerations as a function of time were measured in combination with a TDC signal 
from the second cylinder, Fig. 4. 

One can conclude from Fig. 4 that: 

— time signal of vibration has the highest value near the top dead centre, which means that that 

part of the signal is generated by the closing of delivery valve (see Fig. 1.), 
— instantaneous values of that signal part are not constant, they change from cycle to cycle. 


Fig. 2. Delivery and suction valves [2]: 1 — seat, 2 — plate, 3 — plate spring, 4 — damping plate, 5 — stop plate, 6 — 
guiding ring, 7 — bolt 


Fig. 3. A section across compressor components essential for vibration propagation (arrows show points where 
vibration transmitters were mounted): 1 — axes of valves, 2 — cylinder head, 3 — plate of the valve seat, 4 — cylinder 
block, 5 — axis of the piston movement, 6 — seal 
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Certain valve conditions were simulated in order to determine the influence of valve condition 
on the vibration signal. Then vibrations of the delivery valve bolt were measured. Tests were made 
without an external load of the compressor, i.e. with the cylinder cover removed and with low 
rotary speed. Simulated compressor conditions are given in Table 1. The peak-to-peak value of 
vibration was measured in a period much longer than one revolution. An example diagram is 
shown in Fig. 5. 

The peak-to-peak values changing in time lead to a conclusion that after a time needed to 
stabilize the working conditions, the peak-to-peak value for simulated states reaches an average 
(constant value) and variable value A (difference between maximum and minimum values). Both 
values depend on the simulated condition of the valve. Table 2 contains simulated conditions and 
corresponding peak-to-peak values of vibration accelerations of the valve bolt. 

Presented below in Table 3 are the peak-to-peak values of vibrations measured on the valve, 


valve plate and the cylinder head for one “normal” condition of the valve. 
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Fig. 4. Time / vibration acceleration diagrams and a corresponding TDC signal from the second cylinder, signal 
measured on: 1. valve bolt at 72 rpm and head removed, 2. valve plate at 856 rpm, 3. head at 856 rpm 
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Tab. 1. Mea 


surements table 


Item Condition Simulation 
1 | Normal Valve is operational 
2 | Fracture A of valve plate Radial "natural" fracture along 1/3 of the radius 
3 | Fracture B of valve plate Lack of 1/3 external ring of valve plate 
4 | Fracture C of valve plate Radial fracture along the whole radius length 
5 |High spring stiffness 3 springs 
6 |Low spring stiffness 1 spring 
7 | Suspension of valve plate No valve plate 
8 |Spring wear No spring 
9 |Backlash of two connected components | Loosened valve nut 
Backlash in valve-plate connection Loosened tightening ring 


TL AUNT T 


Fig. 5. Peak-to-peak values of vibration accelerations of an operational valve 


Tab. 2. Measures of the time function of peak-to-peak values of vibration accelerations (valve bolt) 


Condition Constant value (mean) [m/s?] А [5/2] 

Normal 79.2 38.7 
Fracture A 61.6 37.5 
Fracture B 4.15 (uncertain value) 
Fracture C 55.6 34.4 
High stiffness of spring 71.4 25.0 
Low stiffness of spring 44 3.0 
Suspension of valve plate 7.55 6.6 
Spring wear 324 45.2 
Backlash of two connected components 39.35 15.55 
Backlash in valve-plate connection 69.975 25.1 
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Tab. 3. Peak-to-peak values of vibration accelerations generated by a valve working in the "normal" condition as a 
time function, for three measurement points 


Measurement point: Constant value (mean) [m/s?] A [52/8] 
Valve bolt 79.2 38.7 
Valve plate 11.4 10.0 
Cylinder head 34 5.1 


5. Conclusion 


Simulated conditions of the valve tested consist in lower stiffness of the spring and lower 
tightness between the plate and t seat. Only a spring stiffness increase causes the valve lift to 
decrease. Lower spring stiffness leads to a drop in the mean peak-to-peak value. In the examined 
case — spring stiffness increase and related valve lift decrease — the decrease of mean peak-to-peak 
value turned out to be larger than the simultaneous increase caused by increased spring stiffness. 
Fractures in the valve plate cause a decrease in the mean peak-to-peak value. The causes may be as 
follows: decrease of the difference value of pressures on both sides of the valve and an increase of 
impact duration (due to lower plate stiffness). It seems that practically the wear of plate valves 
resulting in an essential increase of valve lift is unlikely. Therefore, it seems impossible that the 
vibration signal value will increase due to valve wear. In case of plate valves, one should expect 
that values of vibration measures will decrease depending on its technical condition. Additionally, 
these values depend on the measurement point: the signal weakens as the distance to a measuring 
point increases and if there are more indirect components in between. 
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Abstract 


The article provides a discussion on the studies and analysis of exposure to vibrations generated by the motor 
vehicle. For the driving safety and comfort it is very important what kind and values of vibration are to car body. As 
the vibration source the motor-engine was chosen. The experiments were conducted on the car vehicle which was 
placed on the special test racks. It allows eliminate the road roughens impact on the suspension and in result to car 
body. The changes of the vibration signals from the motor-engine, floor panel and seat were observed and measured 
in 3 axes. These vibrations are producing a level of discomfort for driver.. 


Keywords: engine vibration, exposure to vibration, vibration propagation 
1. Introduction 


The vibrations are inseparable phenomena during operating and driving by all means of 
transport. Vehicles in motion are forced to vibration mainly by the road [20]. There are many of 
different vibration sources in vehicles as well. Some of them are main source of vibration during 
operating but without any movement. The motor engine should be considered as one of the most 
important vibration generator when the car doesn’t drive. For the driving safety and comfort it is 
very important what kind and values of vibration are transfer to car body and driver. To improve 
the vibration isolation transferred to car body the identification of all sources of vibration and 
vibration material propagation should be taken into consideration. Thirty ears ago vehicles 
constructions were very heavy but than some policy started to demand to make the fuel 
consumption reduce. This started the evolution in material used in vehicles construction. Currently 
the car industry faces a crucial weight problem resulting from increasing customer demands in 
terms of safety and performance. Car bodies contribute 25% to the total weight of a car. Light 
weight metals are solution to decrease the body in weight. An increasing use of metals such as 
aluminium and magnesium in the automotive industry shows that there is still large scope for 
improvements [4,16,18]. It is very important to test new materials in many aspects. One of the 
more important is to guarantee the founded vibration damping level in terms of safety and comfort. 

This paper take into account exposure to vibrations generated by the motor vehicle as the 
detrimental effect on the safety and comfort in mean of transport. The human response to vibration 
is depending on the values, frequencies and directions. The driver and passengers exposure to 
whole-body vibration of the vehicle can affect from short-term body discomfort and inefficient 
performance to longterm physiological damage. 
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2. Human perception of vibration in vehicles 


Motor engine should be considered as the vibration generator as well. This kind of machine 
generate a disturbing force of one sort or another, but the frequency of the disturbing force should 
not be at, or near, a natural frequency of the structure otherwise resonance will occur, with the 
resulting high amplitudes of vibration and dynamic stresses, and noise and fatigue problems. 
Rotating machinery such as motors can generate disturbing forces at several different frequencies 
such as the rotating speed and blade passing frequency. Reciprocating machinery such as 
compressors and engines can rarely be perfectly balanced, and an exciting force is produced at the 
rotating speed and at harmonics. There are two basic types of structural vibration: steady-state 
vibration caused by continually running machines such as engines, air-conditioning plants and 
generators either within the structure or situated in a neighbouring structure, and transient 
vibration caused by a short-duration disturbance such as a lorry or train passing over an expansion 
joint in a road or over a bridge. 

Ride comfort is extremely difficult to determine because of the variations in individual 
sensitivity to vibration. There are some research result published on effects of whole-body 
vibration and the ride comfort limits [6,9,10]. Some studies on heavy vehicles ride comfort mainly 
focus on road vehicles running on the ground [8,11]. There are some papers describing how 
vibration interferes with people's working efficiency, safety and health [1,14,15]. Therefore many 
researchers have concentrated their efforts on reducing the amount of vibration from vehicles. 
There are many reports describing the measurement of the transmissibility of the human body 
under vibration [7,12,13,17]. Some interesting researches were conducted for the low frequency 
discomfort for human analysed [19]. Ride vibrations are transmitted to the driver buttocks and 
back by the seat. The floor panel, pedal and steering wheel transmit additional vibrations to the 
feet and hands of the driver. These vibrations are producing a level of discomfort for driver. 

Human perception of vibration is very good. It is a real challenge in structural design to ensure 
that the perception threshold level is not exceeded. Polish PN-91/N-01354 standard specifies 
methods for assessing exposure to vibrations of the overall impact on the human body. The 
parameter value can be used vibration dominant effective weighted vibration acceleration 
determined from the formula: 

Aj, = MaX|L.4- RMS(a,),1.4- RMS(a ,),RMS(a,) }, (1) 
where: 
Awmax - Vibration dominant effective weighted vibration acceleration, 
a, — acceleration of vibration in X axis, 
a, — acceleration of vibration in Y axis, 


а. — acceleration of vibration in Z axis. 
3. Research 


Under the studies in question, active experiments were undertaken featuring measurements of 
vibration accelerations in a three directions in three selected points to analyse propagation of 
vibration generated by engine to driver feet and back.. It were recorded the vibration in three 
orthogonal axes (X,Y,Z). The purpose of the research was analysis of the car body vibration 
generated from motor engine. The experiments were conducted on the car vehicle which was 
placed on the special test racks. It allows eliminate the road roughens impact on the suspension 
and in result to car body. The paper presents some results of measurements vibration of motor 
engine, floor of the car under the driver foots and driver seat. It enables to analyse the way of 
vibration transfer from the source to driver. Ride vibrations are transmitted to the driver buttocks 
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and back by the seat. The floor panel transmit additional vibrations to the feet of the driver. These 
vibrations are producing a level of discomfort for driver. 


Fig. 1. Measurement of vibration of combustion engine in 3 directions 


The established scope of research enables to observe changes of the vibration for chosen points 
on the vibration propagation way from engine. The three orthogonal axes were analysed 
separately. The comparison of the acceleration of vibration signals allows determine which 
directions of the vibration propagation is parent. The proposed methodology allows estimate 
influence of vibration generated by engine on human perception of vibration. The charts in this 
section illustrate time realization and spectrums of the vibration signals recorded during the 
experiments. The figures are grouped by the measurement points and contain vibrations in 3 axes. 


Vibration signals - X 
- neutral gear - engine block, X axis 


o 


acceleration [m/s?] 


B 


40 L П П П П П П 1 П 


time [s] 


Vibration signals - Y 
- neutral gear - engine block, Y axis 


acceleration [m/s?] 


time [s] 


Vibration signals - Z 
- neutral gear - engine block, Z axis 


acceleration [m/s?] 


time [s] 


Fig. 2. Engine block vibration signal recorded in 3 axes 
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Spectrum of the vibration signals - X 
- neutral gear - engine block, X axis 
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Fig. 4. Floor panel vibration signal recorded in 3 axes 
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Spectrum of the vibration signals - X 
- neutral gear - floor panel under the driver foot, X axis 
3000 rpm 
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Fig. 5. Spectrums of the floor panel vibration signal recorded in 3 axes 
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Fig. 6. Driver seat vibration signal recorded in 3 axes 
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Spectrum of the vibration signals - X 
= neutral gear - driver seat, X axis 
3000 rpm 
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Fig. 7. Spectrums of the driver seat vibration signal recorded in 3 axes 


To compare the energy of the vibration signal in orthogonal axes and energy loss during 
propagation of vibration the root mean square (RMS) was calculated. The results of RMS 


calculation have been depicted in Fig. 8. 
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Fig. 8. The RMS values of the vibration in measured points and in 3 axes 
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This simple signal processing allows calculation of the vibration dominant effective weighted 
vibration acceleration. This estimator is used for assessing exposure to vibrations of the overall 
impact on the human body. The results have been depicted in Fig. 9. 
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Fig. 9. А „ах vibration dominant effective weighted vibration acceleration for the engine vibration propagation 


4. Conclusion 


The results of the conducted research allow detail analysing exposure to vibrations generated 
by the motor vehicle. As it has been presented in Fig. 2 and 3 the vibration of engine block are 
quasi stationary. The largest acceleration values can be observed for the longitudinal vibration (X 
axis). It can be effected by the engine mounting elements, which characteristics are mostly 
designed for vertical displacement preventing. Signals registered on floor panel and the seat have 
different dynamics structures. It can be observed domination of lower frequencies in the signal. 
Based on the results of the vibration dominant effective weighted vibration acceleration it can be 
assumed the vibration generated by the engine are strongly absorbed during the structural 
propagation and the human exposure to vibrations has small values. 
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Abstract 


Paper presents results of an experimental verification of the method of quantitative evaluation of engine operation, 
presented in the literature, exemplified by a low-power internal combustion piston engine. In accordance with such 
interpretation, engine operation may be presented as a physical quantity defined as operation indicator. The paper 
presents results of preliminary tests, processed in that aspect, carried out on an engine test bench. The results have 
been used as a classifying measure of the engine reliability state. 


Key words: operation, internal combustion piston engine, reliability 
1. Introduction 


Rational operation of technical devices and systems requires right operational decisions to be 
taken. The decisions are based on various premises and of a significant importance among them 
are reliability characteristics of the operated object. 

The reliability analysis is a complex procedure of processing the empirical results gathered 
in the operational practice supplemented with analytical considerations, which procedure requires, 
apart from other elements, having access to comprehensive information on the current technical 
and power state of the object in question. 

For identifying and predicting the technical states a suitable diagnostic system must be used. 
Without such a system it is impossible to acquire the information necessary for controlling the 
operation process in general and predicting the reliability level in particular. 

Due to obvious technical and economic reasons, not all the possible technical states can be 
diagnosed even in the case of very simple technical objects. 

As the process of changes of technical state is stochastic, continuous in terms of the states 
and time, it is necessary to subdivide that infinite set of states into a finite number of subsets 
(classes), clearly and permanently identifiable by means of the functioning diagnostic system. In 
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terms of reliability, selecting from the set of all possible technical states the subsets of significant 

states is connected with the need of defining precisely the broadly understood failure. For the task 

reliability (associated with the performed task) of complex technical objects, as are the devices and 

functional subsystems of a marine power plant, it is useful and practical to distinguish at least two 

groups of failures understood as: 

e events causing inadmissible deviations of the system (device) operation parameters from the 
critical values — major failures causing unfitness for performing any task whatever; 

e events causing only deterioration of the operating qualities — minor failures causing in effect 
higher cost of performing the task or a need of modifying it. 

By adopting such subdivision criteria — capability of an object of carrying out the required 
functions connected with the objectives defined by the operating system in given conditions — a 
simple and practically useful for the operation process control set S of classes — technical states 
can be distinguished: 

a) State s; — subset of technical states of full task fitness. An object in that state 1s 
capable of performing all the tasks it has been prepared for in the design and 
manufacturing phase and the values of its operation (power, economic, ergonomic 
etc.) parameters are kept within the admissible range. 

b) State s; — state of partial unfitness (partial fitness) for task execution. Including e.g. 
an engine to that class of states should occur in the case of: 

e incapability of fulfilling all the required ship design parameters, e.g. the specified 
speed or sailing range, 

e incapability of maintaining the manufacturer-guaranteed power and economic 
indicators of the operational process effectiveness, e.g. unit fuel consumption. 

c) State s; — state of full task incapability which makes it impossible to use the object in 
accordance with its intended purpose. This subset includes all the states characterised 
by very significant destruction (extensive failure) of the object. 

In the case of transport power systems, the internal combustion piston engines take particular 
place. The universality of their use, impact on the system safety and ecological properties and also 
the connected operational costs make the knowledge of their current technical state particularly 
important from the point of view of timely reaction to any irregularities. 

Among many indicators allowing to include that technical state in a synthetic way to one of 
the above defined classes, it seems sensible to consider the engine operation in such a valuating 
manner that the operation can be determined simultaneously by energy and time. 

Operation within a [0, t] time range can then be interpreted as a physical quantity — hereinafter 
called operation indicator (D) — determined by a product of the energy E = f(t) variable in time and 
the time, which may be generally expressed by the following relation [1, 4]: 


p= [£y (1) 


In the case of general analysis of the compression-ignition engine operation it may be assumed 
that the energy generated by combustion of fuel in the cylinders produces the engine torque. An 
effect of transmission of the torque from engine to receiver is the work Le, which in this case may 
be determined from the expression: 


Le=M,-2an-t (2) 
where: 
L, — useful work, 
M, — mean torque, 
n — engine rotational speed. 
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From relations (1) and (2), the engine operation indicator may be determined by the formula: 
t 
D -2z M,ntdt (3) 
[U 


By introducing to further considerations the following notions: 

e required operation — Dy necessary for carrying out the intended task (e.g. a transportation task 
— sea transport of cargo within a determined time — which means maintaining the determined 
average ship speed, i.e. average power developed by the ship main propulsion engine(s)), 

e possible operation — Dy which the engine in a given technical state and in given functioning 
conditions is capable of performing, 

and by verifying the relation [4]: 

Dyu? Dw (4) 

a criterion can be obtained of evaluation of the fitness for use, e.g. in accordance with principles 

presented in [4]. 

Further in this paper an experimental verification is presented of the above mentioned 
possibility of using the operation indicator as a general classifier of engine state, based on the test 
bench investigation of a low-power compression-ignition Robin-Subaru EX17 engine. 


2. Engine diagnostic tests 


The investigation of a Robin-Subaru EX17 engine was carried out on a laboratory test bench 
(Fig.1). This is a single-cylinder four-stroke carburettor engine fuelled with lead-free petrol, 
splash-lubricated with the SAE 30 CastelGarden oil. 


Table 1. Basic characteristics of the Robin-Subaru EX17 engine 


Rated power 2.6 kW 
Cylinder diameter 67 mm 
Nominal rotational speed 3000 revs/min 
Engine displacement 169 cm? 
Compression ratio 85 


The engine power was determined by means of an electric rotational brake with a control 
system allowing to measure the engine torque M [Nm]. 

The rotational speed was measured with an inductive gauge cooperating with a toothed 
wheel, where the consecutive teeth are the shaft position indicators. 


a) 


Fig. 1. View of the Robin-Subaru EXI7 engine test stand (a) and the engine cross-section (b) 
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The investigation was performed in the form of a simple active experiment and was aimed at 
verification of the operation time of engine in three different fitness states, operating on 100 cm? of 
the 95 Lotos lead-free petrol with 720 kg/m? density and with the same loads. The following state 
classes were distinguished: 

State | — state of full fitness. 

State 2 — state of partial fitness — damaged lubrication oil splasher (Fig. 2). 

State 3 — state of partial fitness — dirty air filter (Fig. 2). 


Air filter 


Fig. 2. Damaged oil splasher (a) and location of the air filter (b) 


3. Investigation results and their analysis 


The test was carried out for three values of the engine braking torque: 2.5 Nm, 5 Nm and 
7.5 Nm, the engine operation time top on 100 cm? of fuel was measured for each of those loads. 
Each measurement was repeated five times. Examples of the measurement results are presented in 
Fig. 3 and 4. 
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Fig. 3. Engine running time with braking torque value 2.5Nm 


34 


Torque [Nm] 


engine technically proficient 
А . 


А А 


Тіте [5] 
Fig. 4. Engine running time with braking torque value 7.5Nm 
The analysed engine operating conditions made it possible to determine: 


e numbers of the performed engine operation cycles within the top time for different states of 
fitness and different loads [2,3] (Table 1), 


Table 1: Comparison of the numbers of engine operation cycles 


Number of engine Number of engine Number of engine 
Engine state cycles with load cycles with load cycles with load 
M-2.5Nm M-5Nm M=7.5Nm 


State 1 323434 366504 297609 
State 2 294742 296688 232500 
State 3 172500 227050 207000 


e values of work performed by the engine in the period (0, top) of its operation in accordance 
with expression (2), 


Load - 7 Nm 
400000 — — Engine state 1 24] 
2 
Engine state 2 Pe 
gm - 
350000 = - -Engine state 3 Él »- 


Work in the period (0, t) [1] 


Time [s] 


Fig. 5. Values of work performed by the engine in the period (0, top) of its operation with a 7.5 Nm torque load 
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Fig. 6. Values of work performed by the engine in the period (0, top) of its operation with a 5 Nm torque load 
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Fig. 7. Values of work performed by the engine in the period (0, top) of its operation with a 2.5 Nm torque load 


e values of the engine operation indicators in the period (0, top) of its operation in accordance 
with expression (3). 
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Fig. 8. Values of the engine operation indicator in the period (0, top) of its operation with a 7.5 Nm torque load 
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Fig. 9. Values of the engine operation indicator in the period (0, top) of its operation with a 5 Nm torque load 
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Fig. 10. Values of the engine operation indicator in the period (0, top) of its operation with a 2.5 Nm torque load 


Analysis of the results presented in figures 5 to 10 and in Table 1 allows to formulate the 
following conclusions: 

a) engine operating in the state of full fitness is characterised by the longest operation time and 
the greatest number of performed work cycles with a given braking torque; 

b) assessment of the engine fitness by the value of work performed in the period of its operation 
(figures 5 to 7) may lead to erroneous conclusions, because, as it can be seen in Fig. 7, work 
performed by the engine in a partial fitness state (їп this case state 2) reached a higher value 
than in the full fitness state (state 1); 

c) using the engine operation indicator as a classifier of the engine reliability state always leads to 
correct conclusions, because for the full fitness state it takes the highest values regardless of 
the load. 


4. Final remarks and conclusions 


Presented interpretation of engine operation makes it possible to formulate a preliminary 
assessment of the engine degree of fitness. 

The presented method is, in the opinion of the Authors, a valuable complement to the so far 
applied methods of describing the reliability aspects of a so important element of most drive 
systems as an internal combustion piston engine. The basic advantageous feature of the method is 
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connecting the evaluation of power with the time of task execution, which is very important 
particularly in the case of carrying out long-lasting transportation tasks. 

An additional advantage is its universality and therefore applicability to the reliability 
analysis of any power device or subsystem. 
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Abstract 


The paper presents applicability of the theory of semi-Markov processes to determine a limiting 
distribution of the process of changes of technical states of fuel systems for marine engines running on heavy 
fuel oils. The proposed study of this process includes the components of such fuel systems like: 1 - injectors, 
2 - high pressure hoses, 3 - injection pumps, 4 - low pressure hoses, 5 — fine filters, 6 - coarse filters, 7 — 
fuel-feed pump, 8 — fuel heater and 9 — viscosistat with a viscometer. A semi-Markov state transition model 
consisting of ten states has been developed for such systems. Application of technical diagnostics has been 
pointed to be necessary to investigate the state transition process for the systems. The conclusions presented 
furthermore in the paper provide advantages which (according to the author) are of the most importance in 
the design and operation phases of fuel systems in marine diesel engines that run on heavy fuel oils. 


1. Introduction 


During operation of internal combustion piston engines, the fuel combustion process 
which proceeds in marine diesel engines is particularly important. This process should run in 
a way that provides as low fuel consumption as possible. It is difficult to reach such a 
combustion process in internal combustion engines running on heavy fuel oils. This is due to 
the need of heating the fuel to a high temperature to be able to produce in the engine working 
spaces a fuel-air mixture having a proper micro- and macro-structure. Heating the fuel is 
indispensable to reduce the viscosity of heavy fuel oils which should range from 15 + 25cSt 
(for low speed engines), and 12+16cSt (for medium speed engines) [11, 12]. Such a structure 
can be gained when the engine power supply system (engine fuel system) for heavy fuel oil 
finds itself in a state of full ability (so). This means that all components of the system must 
stay in ability states. Such components include [11, 12, 14]: 1 - injectors, 2 - high pressure 
hoses 3 - injection pump, 4 - low pressure hoses, 5 - fine filters, 6 - coarse filters, 7 — fuel-feed 
pump, 8 — fuel heater and 9 — viscosistat with a viscometer [11, 12]. Failure of any of the 
components in the engine fuel system makes the component disable. This causes a failure and 
inability state of the entire fuel system. Therefore, the users of marine heavy fuel engines are 
interested in what the probability is that at any moment in a long time of operation 
(theoretically t — oo) the engine fuel system is in the state of ability (so) or states of inability 


(50 : 51, $5, ety So). 39 


2. General identification of state changes for marine diesel engines 


The state of ability (so) of each fuel system in marine internal combustion engine is 
found (as mentioned in the introduction) when all components in the system stay in ability 
states. If any of the components gets a failure, the fuel system must be considered to be 
damaged and, therefore, in state of inability (5, ). A failure of the i-th (any) of the components 
in this system is fixed by the engine user during performance of a maintenance service. As a 
result of this action the faulty component is renovated and in consequence stays no longer in 
the inability state s; € Ld (i = 1, 2, 3,...., 9), instead of this, it takes a condition that provides 
an ability state so to the entire fuel system. Thus, a marine engine fuel system in the phase of 
operation can take one of the states from the set: 


S = {50, 51, $2, ..., S9} (1) 


with the following interpretation: 

e 50 — ability state of a fuel system, 

e (si i = 1, 2, 3,...., 9) — inability state of a system due to a failure of, respectively: 
injectors (si), high pressure hoses (52), injection pumps (5з), low pressure hoses (54), 
fine filters (sc), coarse filter (sc), fuel-feed pumps (s7), fuel heater (ss) and viscosistat 
including a viscometer (so), i.e., that ses 

Therefore, the process of occurring states sie S (i = 0, 1, 2, 3,...., 9) in sequence can be 
observed. The operational practice shows that the unconditional duration 7; of any state s;, and 
Ту of the state s;, provided that the state s; is the next one, are random variables [2, 3, 7, 8]. 
Furthermore, it is obvious that the system can stay in these states with a determined 
probability P;, while a transition from the state s; to s; proceeds with the probability ру [3, 8, 
10, 13]. 

It can be assumed that a fuel system in a new engine is in the ability state (so). The same 
state of it can be found after its renovation made during performance of an appropriate 
maintenance service [3, 7, 12]. Due to the fact that the time of renovation performance is a 
random variable, the occurrence of the state so 1s a random event. The system may stay in this 
state for the time period To with the probability Po. The time period To is a random variable, 
because the value of the time results from the occurrence of a failure in the system, which is a 
random event. A state of inability of the fuel system occurs as a result of its failure. This state 
is a consequence of occurrence of disability of any component within the system. Thus, for 
this system there can be distinguished the following inability states: s;e S (i = 1, 2, 3,...., 9). 
When the cause of the inability is a failure of injectors, the state transition from so to s; takes 
place after the time period То with the probability poi. The occurrence of the state s; will 
cause a need to implement some service in order to renovate the system and regain its state so. 
The fuel system returns to the state soy after the time period То with the probability pio. 
However, if the cause of the state of inability of the fuel system 15 a failure of high pressure 
hoses, the transition proceeds from the state so to 52. This transition takes place after the time 
period 702 with the probability po». In this case, the system can return to the state so after the 
time period 75o with the probability p20, provided that a proper service is performed. The same 
situation applies when the fuel system loses so in the event of a failure of any of its other 
components. For each of the considered cases the conclusion can be made that the current 
state of a fuel system and the time period of its duration depend significantly (even 
exclusively) on the immediate previous state. This is obvious when the consideration includes 
the fuel system states occurring in sequence: s;e S (i = 0, 1, 2, 3,...., 9). For example, if a fuel 
system is in the state so, with the time of its failure this state transits to one of the possible 
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states of sie S (i = 1, 2, 3,...., 9), depending on which of the components of the system is 
damaged. When the inability state of the fuel system is 51, then after its renovation and return 
to the state so, just this state (so) and the time period of its duration will depend only on the 
state 51, not on those states that were earlier, so the states of sie S (i = 2, 3,...., 9). Duration of 
the state so, after implemented service will in turn depend on how well the service restoring 
this state (so) was performed and also on how much the system will be loaded during 
operation [5, 6, 7, 12, 14]. The above considerations show that a semi-Markov model of a real 
process of changes of its states can be considered in case of fuel systems for marine diesel 
engines running on heavy fuel oils. Such a model is possible to be built because the process 
states can be determined in such a way that the duration of the state existing at the instant т, 
and the state obtainable at the instant т„ do not depend stochastically on the states that 
occurred earlier or duration of these states. The theory of semi-Markov processes can be 
applied to build a semi-Markov model (W(t): t 2 0} of the real process of changes of technical 
states for the operation phase of propulsion systems. The characteristics for the model are [7, 
8, 10, 12, 13]: 

1) the Markov condition is satisfied so that in future the evolution of the process of 
changes of states of a fuel system during its operation phase, for which the semi- 
Markov model has been built, would depend only on its state at the given time, not on 
the system operation in the past, thus so that the future of this system would not 
depend on the past but on the present, 

2) random variables 7; (denoting the duration of the state 5, independently which state it 
will be followed by) and Ту (denoting the duration of the state ,,s7’, on the condition 
that the state "зу" will be the next one in the process) have other than exponential 
distributions. 

Modeling, which allows to develop a semi-Markov model of the process of changes of 
fuel system states includes the outputs of the conducted analysis of the changes in these states 
[2, 5, 8, 12, 14]. The states have been recognized to be the most important in studying the real 
process of changes of technical states of fuel systems for marine diesel engines on heavy fuel 
oils. 

The presented description of the changes of technical states for fuel systems is a result 
of the noticed fact that for these systems, just like for other functional systems in diesel 
engines, prediction of changes in their technical states can be made when the current states of 
the systems and the conditions under which they are expected to operate in future are known. 
This also applies to many other technical objects, for example such as gas turbine engines, 
reciprocating compressors, centrifugal pumps, etc. [2, 4, 7, 12, 14]. This fact which also 
indicates that changes of technical states of this kind of objects are not closely related to the 
time of their operation (work), is explained as a scientific hypothesis by the publications [2, 
7]. With regard to the fuel systems investigated herein, this can be explained in a form of the 
following scientific hypothesis (H): the process of changes of technical states of any fuel 
system (defined as a random function whose the argument is time and the values are random 
variables denoting the current technical states), running under rational operation (that is 
operation based on economic calculation) is a process with asymptotically independent 
values, because any particular state of the process investigated at any instant tn (n = 0, 1, 
es M; TO € U <... € Tm), depends significantly on the immediate previous state, not on the 
states that occurred earlier or the periods of their duration. 

It should be noticed that the formulated hypothesis does not contain any such 

contradictions that might falsify it in the logical sense, still before testing. The consequences 
of this hypothesis are as follows [2, 7, 8]: 
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1) Ki — probabilities (pi; i +j; ij Е N) of the process transition, i.e. changes of fuel 
system states from any state ,,s;”, which the process (system) is in, to the next (any) 
state " s;", do not depend on the states which the process was earlier in, 

2) K2— periods of unconditional duration of the particular states "s;" of the process of 
changes of fuel system states are stochastically independent random variables (T; i € 

N), 

3) Кз — periods of duration of each of the possible-to-occur states "s;" of the process of 

changes of fuel system states, provided that the next state is one of the remaining 
process states "s;", are stochastically independent random variables (Ту; i +j; ij Е N). 
The presented consequences reveal the probabilistic law for the changes of technical 
states of any fuel system. They are not mutually contradictory, and their logic truth raises no 
doubts. Thus, the non-contradiction condition for the consequences is satisfied. This means 
that nothing stands in a way to make advantage of the mentioned consequences regarded 
jointly as one consequence (K) in order to test the hypothesis (Н) empirically, so to make its 
verification if should be accepted or falsified. This verification consists in experimental 
testing of the truthfulness of the consequences К; e К, i = 1, 2, 3. А non deductive inference 


method called a reductive inference can be applied to verify this hypothesis. 


3. Semi-Markov model of the process of changes of technical states of fuel 
systems for marine diesel engines 


Application of a semi-Markov model of the process of changes of technical states of a 
fuel system for a heavy fuel engine allows to: 

* take into account preventive maintenance required to restore the system which 
becomes completed upon completion of maintenance of particular system's 
components and 

e investigate more than two reliability states of an engine fuel system as well as its 
components. 

This allows to investigate a semi-Markov model of changes of technical states of fuel 
systems for marine internal combustion engines with a ten-element set S (1). This means that 
the limiting distribution of the process of changes of ability and inability states of heavy fuel 
supply systems for marine diesel engines, can be defined by using the model of the process of 
changes of the system states, developed in а form of a semi-Markov process ((W(t): t > 0} 
with a set of states 5 = s; i = 0, 1, 2, ... 9. Interpretation of the states s; є S(i = 0, 1, 2, ... „9) 
in accordance with the description in equation (1) is as follows: so— ability state of a fuel 
system, s; — inability state of a fuel system due to a failure of an injector (injectors), 52 — 
inability state of a fuel system due a failure of a high pressure hose (hoses), s3 — inability state 
of a fuel system due to a failure of an injection pump (injection pumps), s4— inability state of 
a fuel system due to a failure of a low pressure hose (hoses), 55 — inability state of a fuel 
system due to a failure of fine fuel filters, 56 — inability state of a fuel system due to a failure 
of coarse filters, s; — inability state of a fuel system due to a failure of a fuel-feed pump (fuel- 
feed pumps), ss — inability state of a fuel system due to a failure of a fuel heater, so — inability 
state of a fuel system due to a failure of viscosistat or viscometer. 

Graph of state transitions for the process (W(t): t2 0} is shown in Fig. 1 

Transitions of the states s; (i = 0, 1, 2, ..., 9) occur in succeeding times /, (n € N),, where 
at the time to = 0, the engine fuel system is in state sọ. The state s, lasts until a failure of any of 


the components of the engine fuel system. However, the states (si(i = 1, 2, ..., 9) last until the 
damaged engine fuel system is restored. 
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The above considerations show that it can be assumed that the state of an engine fuel 
system at time ¢,+; and the time period of duration of the state at a time ¢, do not depend оп 
states that occurred at times fo, й, ..., fn-1 or time periods of their durations. Therefore, the 
process { W(t): t > 0} is a discrete-state continuous-time semi-Markov process [2, 3, 4, 8, 10, 


13]. 
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Fig. 1. Graph of state transitions for the process (W(t): t 2 0}: 50 — ability state of an 
engine fuel system, {so} = set of inability states of an engine fuel system: 


{50} = %&],82,83,.589 1, si € S(i = 1, 2, ... ,9) — states with the interpretations respectively: s] — 
inability state of an engine fuel system due to a failure of an injector (injectors), s2 — inability state of an engine 
fuel system due a failure of a high pressure hose (hoses), 53 — inability state of an engine fuel system due to a 
failure of an injection pump (injection pumps), 54 — inability state of an engine fuel system due to a failure of 
a low pressure hose (hoses), s5 — inability state of an engine fuel system due to a failure of fine fuel filters, 
56 — inability state of an engine fuel system due to a failure of coarse filters, s7 — inability state of an 
engine fuel system due to a failure of a fuel-feed pump (fuel-feed pumps), s8 — inability state of an engine fuel 
system due to a failure of a fuel heater, s9 — inability state of an engine fuel system due to a failure of viscosistat 
or viscometer, Tij - time period of the process duration in the state si provided that the state sj is the next 
one, pij = probability of staying the process in the state si, provided that the state sj (i, j = 0, 1, 2, ..., 9; iz j) is 
the next one. 


An exemplary realization of the process (W(t): t > 0}, depicting occurrence of technical 
states of a fuel system in any heavy fuel engine when operated, is shown in Figure 2. 

Obtaining (obviously in approximation) the values of the probabilities Pj (j = 0, 1, 2, 3) 
that make a limiting distribution of the process of changes of ability and inability states of a 
heavy fuel supply system for a marine diesel engine requires determination of the initial 
distribution of the process (W(f): t > O; and its matrix function. 

The initial distribution of the process is: 


r йа i=0 
PWO) = 5,} =} | (2) 
\0 dla i-1, 2, ..., 9 


and the matrix function has the form: 
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Fig. 2. An exemplary realization of the process (W(t): t 2 0] for the investigated fuel system components: s0 
- ability state of the components and the system at the same time, si(i = 1, 2, 3, ... ,9) — states as described in 
Figure 1. 


The matrix function Q(# is a model of changes of technical states of a engine fuel 
system. Nonzero elements Qj(t) of the matrix Q(t) depend on the distributions of random 
variables which are the time intervals of staying the process {W(t): t > 0} in the states s; є S(i 
= 0, 1,...‚ 9). The elements of the matrix function Q( are the probabilities of the process 
transitions from state s;, to state s; (si, s; € S) in time no longer than ¢, which are defined as 
follows [3, 8]: 


Qi(t) = Рт) = sj, Trt — Tn < И) = si) = pF iO) (4) 
where: 
Pj — one-step transition probability of the homogeneous Markov chain; 
whereas py = P(Y(v,«1) = Y(t) = sij=limQ, (1) Н 


F(t) — distribution of ће random variable Те) denoting the duration of the state s; of 
the process { W(t): t > 0}, provided that the state s; is the next one in the process. 
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The matrix P of transition probabilities of the embedded Markov chain to the process, 
as this results from the matrix function Q(t) (2), is composed of the elements which are [2, 7, 
8]: Q(t) = pio = 1(/ = 1, 2, ..., 9) and Qy(t) = ро, because Q(t) is a stochastic matrix. 

The process (W(t): t 2 O}is irreducible [1, 3, 4], and the random variables Тез take 
finite positive expected values. Therefore, its limiting distribution [7, 8] 


P, = lim P,() = lim PAW =s,), s; € SG 70,1, 4) (5) 


is of the following form [8]: 
x ,E(T;) 
pou (6) 
È TECT) 
1-0 


The probabilities ry(j = 0, 1, 2, ..., 9) in the formula (6) are limiting probabilities of the 

embedded Markov chain to the process (W(t): t 2 0}. However, E(7;) and E(7;) are expected 
values of the random variables of 7; and 7; respectively, which are the time periods of staying 
the system in states: s; and s; respectively, regardless of which state will be next. 
To determine the limiting distribution (6) it is required to solve the system of equations that 
comprise the said limiting probabilities ry(j = 0, 1,..., 9) of the embedded Markov chain and 
the matrix P of the probabilities of transition from the state s; to s;. Such a system takes the 
following form: 


[ло,ли,л2,73,74,75,76,77,л8,л9] = [ло TNT, NANS, M6, M8, 9]: 


Ут, =1 (7) 


As a result of solving the system of equations (7), by using the formula (6), the 
following relationships can be obtained: 


E(T, ЕТ ЕСТ. 

p= ( 0) Е Ре) (Т,) ҮР, Ре» (T) , 
E(T,)+ > PonElT,) ЕСТЬ) + ? po ET.) = E(T,) + > ро›Е(Т) 

Ba Pos El) Ba Pos ET.) „B= Pos (Г) Р 
E(T,) + 2 poo E(T) E(T,)+ 2 pan ET) E(T,)+ 2 poo EQ) 

Pes E(T5) p o5 E(T7) 

Р, = | 9 s477 9 , = PE) 

E(To)+ Ż,PovE(T,) E(T,)+ 2 Pen E(T,) ЕСТ.) + > poo EQ) 
ЕСТ. 

P, = Pos ( 9) (8) 

E(T,)+ > Pow E(T) 
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Ре is a limiting probability that the engine fuel system stays in the state so (ability state) 
in a longer period (time interval) of operation (in theory at / — оо). Thus, this probability 
defines a coefficient of technical availability of the system for operation. However, P(j = 1, 
2,..., 9) are limiting probabilities of existence of the states s; e S of the systems at £ — оо, thus 
they are the probabilities that its components (and simultaneously the entire system, due to its 
serial reliability construction) are in states of inability. 

Obtaining the probabilities Pj = 0, 1, 2, ..., 9) which make a limiting distribution of the 
process of changes of ability and inability states of fuel supply systems for marine heavy fuel 
diesel engines requires, as it follows from (6) and (8) formulas, an estimation of the values of 
the probabilities ру and expected values Е(7). Obtainment of the values (obviously, in 
approximation) is expensive and labor consuming. This is due to the fact that collection of 
necessary information to estimate pj; and E(7;) requires long-term empirical research using 
appropriate diagnosing systems (SDG) for engine fuel systems [1, 12, 15, 16, 17]. 

Estimation of the probabilities ру and expected values E(7;) is possible after getting 
realizations w(t) of the process (W(t): t > 0) in sufficiently long time period of investigation, 
thus for t e [0, t], whereas the time of the process investigation is: № >> 0. Then, it is 
possible to determine the numbers of: ny(i.j = 0, 1, 2, 3; i +), transitions of the process(W(t): t 


2 0j from state s; to s; in sufficiently long time and to define the values of the estimator P, of 


the unknown probability ру. The estimator of the transition probability ру with the highest 
reliability is the statistics [7, 8]: 


^ Ni . . 
i;== 7 Lj-0,1 2,3, (9) 
YN; 
j 
whose the value р = "i isan estimation of the unknown transition probability p 3 


ij Yn Р ў 
J 

The w(t) course of the process W(t) сап also give the realizations 4”, m = 1, 2, ..., nj; of the 
random variables 7;. Application of a point estimation allows for ease calculation of E(7;) as 
an arithmetic mean value of 4%”. Acquisition of the necessary information in order to 
calculate the probabilities requires employing adequate diagnosing systems (SDG) for the 
mentioned fuel systems which in this case become diagnosed systems (SDN) [5, 7]. 

The presented description of changes of states of a fuel system in a heavy fuel diesel 
engine can be obviously developed when taking into account the states of partial ability of 
such an engine system. 


4. Final remarks and conclusions 


Application of the aforesaid processes in the practice requires to satisfy the two 
following conditions: 
e to collect the mathematical statistics that are necessary to determine the probabilities 
ру and the expected values Tij, 
* to develop a semi-Markov model of changes of technical states of an engine fuel 
system with a small number of its states and simple (in the mathematical sense) matrix 
function Q(t). 
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Application of a semi-Markov process as a model of changes of the specified states of 
an engine fuel system at a defined time (given moment), instead of a Markov process, results 
from that it should not be expected that the random variable Ту) denoting the time period of 
duration of the state s; on the condition that the state s; is the next one and the random variable 
T; denoting the time period of duration of the state s;(i = 1, 2, 3, ..., 9) of an engine fuel 
system, independently of which state will be next, have arbitrary distributions included in the 
set R+ = [0, +оо). Employing a Markov process for this case would be justified if the 
assumption could be made that the random variables 7, and 7; have exponential 
distributions. 

The presented model can be of considerable practical importance because of ease in 
determining the estimators of the transition probabilities ру which are the elements of the 
matrix P (9) and ease in estimating the expected values Е(Т;). However, the considerations 
should take into account the fact that a point estimation of the expected value E(7;) does not 
enable to define the estimation accuracy. Such accuracy is possible by the interval estimation, 
where a confidence interval [tg tj] with random endpoints is determined that contains with a 
certain probability (confidence level) B, the unknown expected value E(7)). 
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Abstract 


This paper presents a modelling of accumulation and dissipation energy in simple drive system. The simple drive 
system is presented by laboratory stand designed for calculate fatigue life of material for bending and rotation tests 
pieces. The laboratory stand was prepared to diagnosis of fatigue life of material by use energetic methods. This idea 
was submitted as utility model to Polish Patent Office. This laboratory stand gives ability to research a influence of 
change of fatigue life for energetic symptoms (like: temperature of material, vibration, Acoustic Emission, rotary 
speed, electrical energy). This paper presents influence of change of fatigue life for rotary speed. The rotary speed is 
value describing effect of accumulation and dissipation energy in drive system (in material). A difference between 
energy taken from power supply and consumed for giving rotary speed is used for fatigue of material. This simple 
theory is presented and modelling in this paper. The next paper (in this book) describes results of research that could 
verify this theory. 


Keywords: modelling, accumulation and dissipation of energy, mechanical drive systems 


1. Introduction 


The necessary of diagnostic of machines is obvious consequence which is result of their 
maintenance. This concerns especially important and expensive machines. For many engineers and 
scientists a very difficult problem is how to calculate or predict the future life of mechanical 
elements that works by high cyclic load. We know that fatigue life is very complicated problem 
and it is difficult to precise calculation. We know too that the fatigue life is effect of accumulation 
of energy in material. The difficult problem is how to calculate this energy. We don't have to do 
this (at work of machines) but we can observe symptoms that are effect of accumulation energy. 
For example a growth of temperature is a value defining a increase of internal energy in material. 
We know that the growth of temperature depends on kind of material (elastic-plastic or plastic- 
brittle) and this is symptom describing the future fatigue life for this mechanical element. 

The fatigue of material at work of mechanical elements could be described by many symptoms 
and parameters. Their are effect of change of properties of material. The cyclic load makes 
hardening or weakening of material that influences on rigidity. This makes that the bearings could 
be less loaded (the lesser deflection of shaft in place of support) and less energy is wasted there. 
The change of rigidity influences on change of vibration of mechanical drive system. For example 
a stable work after hardening of material (after 20 hours of work) makes stability of vibration 
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amplitude. This symptom is useful to diagnostic of technical state all machines or engines. 

The change of properties of material influences on rotary speed. In used laboratory state the 
change of rotary speed is visible diagnostic symptom describing accumulation of energy in 
material. This accumulation of energy describes fatigue life of steel test piece. The observation of 
multi symptoms by use temperature, deflection, vibration, electrical energy: current and voltage, 
Acoustic Emission and rotary speed gives information about effects of accumulation energy in 
material. This same allows to estimate a fatigue life. 

The multi symptoms research allow to make physical and mathematical model. It is necessary 
to verification results and could allow to predict fatigue life. 


2. The purpose of test stand building 


Every mechanical drive system owns a efficiency. The main stream of energy is passed to any 
working equipments. The part of energy is lost in mechanical elements (bearings, shafts, gear 
wheels, oil) which pass energy from source (electrical motor or engine) to working equipments. 
This small part of lost energy is used to change properties of material. The bearing and tooth of 
gear wheels are cyclic deformed. A oil is plastic deformed under cyclic loading of shaft neck. 
These all cyclic changes make growth of temperature and change material's properties. These 
results are effect of dissipation energy in material and dissipation energy in mechanical system 
(growth of vibrations). 

Any delivered energy makes a change of situation of body (substance). This principle concerns 
to every situation: huge machines that are described like closed system and every elements (inside 
of substances). When it is considered to separate elements the dissipation of energy makes change 
of material's properties. Overflow of limit of this dissipation energy in material makes complete 
wear or crack of material. Present research should show how can describe technical state of 
working elements by use theory about dissipation energy in material. This problem is difficult to 
solve but it will give good information about calculate of reliability of mechanical structures. 

The used laboratory stand is prepared to research fatigue life of bending and rotary elements 
working by cyclic loading. This laboratory stand is patterned on well-known fatigue machines - 
Schenck [4] but now it is used to calculate and diagnose of material by use energetic methods. The 
improvement of laboratory stand and new patterns was described in utility model [the number of 
application W.121416]. 


Measurement of 


rotary speed 
Electrical 


motor 


Measurement 
of electrical 
parameters 


Fig. 1. The laboratory stand 
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Acoustic 1 
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Fig. 2. The laboratory stand — place of sensor’s installation 


The principle of laboratory stand is possibility of observation of energetic symptoms for 
bending of rotary test piece. The main principle of laboratory stand is cause that all stream of 
energy is used in mechanical rotary system. The energy had been taken from electrically source 
could be used for two causes: giving kinetic energy to rotary elements or accumulation energy in 
material (what makes his fatigue). This allows to know better symptoms describing change of 
properties of material at cyclic loaded work. The main formula describing behaviour of drive 
model is relationship: 


бам (0 
dt i 


res 


Where: 


J — moment of inertia [kgm], 
œ- angular velocity [rad/s], 

t — time [s], 

Ма- drive moment [Nm], 

Mres — resistance moment [Nm]. 


The principle 1 describes a behaviour of drive system as change between drive moment My and 
resistance moment М. The drive moment is value of energy delivered to system. The resistance 
moment describes wasting of energy in mechanical system: accumulation energy in material 
(increase of internal energy) and dissipation energy in mechanical system (growth of kinetic 
energy — increase of rotary speed). A difference between Mg and Мк is used to give rotary speed. 
So a value and a change of the rotary speed can describe a part of energy that is not used to change 
material's properties. The resistance moment Myes can be described generally as sum of two 
moments: 


M es = M, + M s, Q) 


where: 


Mr- friction moment (in bearings only — the simple model) [Nm], 
Ми — moment of change of material’s structure [Nm]. 
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The Mf is simple function of one physical quantity - loading (stresses in material) because it 
makes more loading of bearings and more friction and lost energy. The principle between M; and 
loading is linear function and can be simple to modeling. In laboratory stand is not more important 
elements where might be lost of energy. 

The moment of change of material's structure М is difficult problem. This is nonlinear 
function of many physical quantities: 

— temperature; 

— deflection; 

— rotary speed; 

— parameters of Acoustic Emission (energy, amplitude, number of events, duration of signal, 

RMS); 

— vibrations. 
The Mgt is following function: 


М = f(T,D,o, AE,V) (3) 
where: 


Mr- temperature [К], 

D — deflection [mm], 

œ- angular velocity [rad/s], 

AE — parameters of Acoustic Emission, 
V — RMS of vibrations. 


3. The physical model of transformation of energy in mechanical drive system 


Three values are necessary to modeling of behavior of drive system: input, output and 
disruptions. The input energy is measured by using electrical parameters: 
— current (measured on rotor); 

— voltage (measured on rotor); 
— electrical input energy (energy delivered from power supply); 

The drive moment Mg is calculated as electromagnetic moment [5]. This principle does not 
contain a efficiency of rotor's system (decrease a voltage in carbon brush and in rotor's winding) 
but this efficiency is high and might have solid number what helps by modeling. The main 
principle is shown below: 


| 4 
NS: a) 
@ 


where: 


U - voltage (on electrical terminal of rotor) [V], 
I — current (on electrical terminal of rotor) [A], 
@ — angular velocity [rad/s]. 


The output values was specified before Ms, = f(T, D, о, AE, V) and present it is not necessary 
precise model describing М; (moment of change of material’s structure). The future research can 
help to make this precise model. In this paper I can show a influence of fatigue life on change of 
the rotary speed w, the temperature T and the drive moment Mg. The main view of all mechanical 
drive system (laboratory stand) is shown below [3]: 
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Fig. 3. The model of laboratory stand (simple mechanical drive system) [3] 


The model on Fig. 3. shows a flow of energy in mechanical drive system. A increase of load 
makes change of diagnostic (energetic) parameters. A increase of load makes decrease of rotary 
speed and increase of temperature. The physical model is shown on Fig. 4. This presents 
relationships between mechanical elements (bearings, absorbers, test piece, electrical motor) and 
physical values that describe accumulation and dissipation energy in mechanical drive system. The 
main principle of physical model is relationship 1. The other relationships are less taken into 
consideration but that are going to model in relationship 3. The disruptions value is the load. The 
load is made by weight hanging under laboratory stand. This makes constant bending stresses on 
all length of test piece. The electrical motor makes that the bending stresses has changeable 
progress (sinusoidal progress). At normal work of machines a change of load is value which makes 
change of behaviour of system (for example as load is increasing as rotary speed is decreasing). 
The description of laboratory stand was presented in earlier articles of Author [1, 2]. 
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Fig. 4. Physical model of simple mechanical drive system and relationships between physical values [3] 
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4. The example of test's results 


The experiment was made for five test pieces. The conditions of work were constant. The 
material of test pieces used in this experiment is designed for mechanical drive elements (shafts or 
gear wheels). The test pieces were made from one round bar that makes smaller spread of results 
what is obvious for fatigue wear. The parameters of starting work and kind of material are shown 
in Table 1: 


Tab. 1. Parameters of test 


number of outside loading starting material 
test piece | temperature rotary speed 
[K] [MPa] [r/min] C45+C 
1,2,3,4,5 296 297 4000 Rm = 830 MPa 
=0,35Rm (round bar) 


The test’s results corroborate a simple theory described in point 3. The progress and change of 
energy and rotary speed is depended on fatigue life. Results of five test pieces are not 
representative population for one test but are a part of bigger research. This test was planned by 
use earlier test (15 test pieces) [2]. A experiment and the laboratory stand were improved (the 
measured values and period of measurement). The behaviour of material under changeable loading 
is known but in present test are more values describing accumulation and dissipation of energy. 

The results of test show a change of behaviour of drive system depending on fatigue life 
(number of cycles until crack). The results confirm the simple theory described in point 3: as a 
fatigue life is higher as rotary speed R(N) is higher and energy from power supply P(N), drive 
moment M,(N) and temperature T(N) are lesser. 

The behaviour of mechanical drive system and accumulation of energy can describe follow: 
energy delivered to drive system is not waste in material but makes higher rotary speed — the 
energy is not accumulated in material and does not make changes in material (elastic-plastic 
deformations and increase of temperature). The temperature is a value describing the internal 
energy. For cyclic loading the internal energy is describing as a range of hysteresis lop between 
stresses and deformations. This statement was described by Szala and Boronski [6] and Kaleta [7]. 
In this test the internal energy is value describing fatigue life and reliability of mechanical drive 
systems. 

The behaviour of mechanical system and principle | is shown on Fig. 5. The Fig 5 makes some 
questions: 

— why does the rotary speed increase when conditions of work are this same? 
— why does the drive moment Ма decrease? 
— why does the resistance moment М; decrease? 


increase 
1 
а E 
Jat =Ma M 
decrease 


Fig. 5. Behaviour of mechanical drive system 


54 


The answers might be contained in progress of measured values. These characteristics show 
change of rotary speed, power supply, drive moment, temperature and current of rotor as function 
of fatigue life. These all characteristics are parameters describing a change of stream of energy at 
work of mechanical drive system. The main results of test are shown in table 2: 


Tab. 2. Results of test 


Test | number of cycles | temperature at rotary speed power supply drive moment 
piece until crack moment of crack | (ending, stable) | (ending, stable) | (ending, stable) 
[-] [K] [r/min] [W] [Nm] 
1 162x10° 313,3 3740 324 0,247 
2 140x10° 315,0 3400 345 0.258 
3 263x10° 308,8 4030 315 0,243 
4 274x10° 309,4 4700 300 0,230 
5 218x10° 312,1 3850 325 0,270 


The current of rotor has a direct relationship with drive moment. For electrical motor 
(permanent current) the drive moment has follow relationship [5]: 


Му=с-ф:1 (5) 
where: 


с — constant factor depending on construction of electrical motor [-], 
9 — magnetic flux [Wb], 
I — current (on electrical terminal of rotor) [А]. 


This relationship shows that drive moment (indirectly energy delivered to drive system) can be 
described by using the current of rotor only (magnetic flux is depended on induction current and is 
permanent). A calculation of drive moment M, by using principle 5 is difficult but a trend of 
characteristic is kept. 


5. Conclusion 


The results of test confirm theory of accumulation energy in mechanical drive system. The 
energy which is not lost in material is used to make increase of kinetic energy (increase of rotary 
speed). This shows that rotary speed is value which can be used to describing fatigue life and 
reliability of drive systems. 

The figure 7 shows interesting progress of power supply. The energetic theory of fatigue wear 
could be described as a energy delivered to material. The endings of progress of power supply 
show Wóhler's characteristic. 
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Fig. 6. Change of physical values in depending on fatigue life (number of cycles N): a) rotary speed R(N), b) drive 
moment Ма, c) power supply P(N), d) temperature T(N), e) current of rotor I(N) 


The results show following conclusions: 
— the change of rotary speed is depending on fatigue life, 


— the temperature is a value describing internal energy accumulated in material and is depending 
on fatigue life, 


— the drive moment and the power supply are value describing energy delivered to mechanical 
drive system and are depending on fatigue life. 


At this test were used other diagnostic parameters: vibrations and Acoustic Emission (AE). 
These parameters confirm energetic methods using by diagnostics of mechanical elements and 
fatigue wear. 

This test allows to predict a behaviour of diagnostic parameters in the next part of experiment. 


Thanks to this test is making now a method of research of the future experiment and a statistical 
analysis. 
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Fig. 7. Wohler characteristic as function of power supply (delivered energy) 


This test allows to predict a behaviour of diagnostic parameters in the next part of experiment. 
Thanks to this test is making now a method of research of the future experiment and a statistical 
analysis. 

The next tests will be doing for higher and smaller loading that allow to observe the parameters 
of drive system (especially a rotary speed) depending on delivered power and a influence on 
fatigue life of material. This all results will be necessary to elaborate a statistical analysis for a 
whole experiment. This allows elaborate characteristic of behaviour of drive system and make 
mathematical model of accumulation and dissipation of energy in this system. 
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Abstract 


The manuscript presents the model of the exhaust gas flow through the exhaust gas duct of the marine 4-stroke 
engine. Presented model are computational fluid dynamic model based on dimensions and the construction of the real 
exhaust gas duct. The measurement parameters from real object are used to the model validation. The simulation of 
the exhaust gas duct throttling by rotational throttling plate was done. Obtained calculation results allow to 
determination of the exhaust gas mass flow for the simulated flow characteristics. The model of turbulence flow was 
taken into account. The gravity forces and the heat transfer phenomena were neglected. Obtained calculation results 
are qualitatively consistent with results obtained from literature. The analyze of the velocity distribution in the exhaust 
gas duct allows to conclusion that the changes of the angular position of the throttling plate causes significant 
disturbances in the exhaust gas flow. The result of this is the decrease of the exhaust gas flow. Additional purpose of 
the manuscript was approximation of the obtained results of the exhaust gas flow for different angular positions of the 
throttling plate. Obtained polynomial function may be useful tool to modeling the combustion process in the engine 
cylinders for the different flow characteristics of the exhaust gas duct. The calculation results allow to determination 
the mass flow of the exhaust gas with mean error equal 11%. 


Keywords: marine diesel engine, exhaust gas composition, toxic emission, CFD model, exhausts gas duct throttling 


1. Introduction 


4-stroke diesel engines are commonly used sources of the mechanical energy in the land and 
the marine applications. Unfortunately they are the sources of the gaseous toxic compounds also. 
Mentioned compounds are products of the hydrocarbon fuels in air. According to regulations of 
the Annex VI to MARPOL Convention [11] all marine engines with nominal power output above 
130kW must apply the nitric oxides emission limits. Emission must be calculated according to 
regulations from “Technical code on control of emission of nitrogen oxides from marine diesel 
engines" [12] and the ISO-8178 standard regulations [3]. According to European Commission 
regulation [2] the diesel engines operated in the on-road vehicles must apply the emission limits of 
carbon oxides, hydrocarbons and solid particles. 

Delimitation of the toxic gaseous compounds emission requires direct measurements of the 
fractions of the mentioned compounds in the exhaust gas during the engine operation with 
specified loads. Quantity of the emission of the toxic compounds requires the calculation of the 
mass flow of the exhaust gas. Three methods are acceptable: 

— direct measurement of the exhaust gas flow, 
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— measurement of the fuel consumption and the air flow to the engine and calculation of 
the exhaust gas flow, 

— measurement of the fuel consumption and calculation of the exhaust gas flow according 
to carbons balance method [3]. 

Direct measurement of the exhaust gas flow in the case of large, marine diesel engines is very 
difficult and not precise. Therefore, often applied method is calculation of the exhaust gas flow by 
carbons balance method. Disadvantage of this method is the calculation error caused by “а priori? 
assumed the chemical composition of the fuel. 

The diesel engines operations causes deterioration of its technical state. Effect of this are the 
changes in the organization of the combustion process in the engine cylinders and changes in the 
emission of the toxic compounds. The solid particles emission causes particles deposition on the 
surfaces of the exhaust gas duct. In [5] author presents the results of the laboratory research about 
influence of the exhaust gas duct throttling on the composition of the exhaust gas. According to 
these results the throttling of the exhaust gas duct causes significant changes of the carbon oxide 
emission. Obtained results allow to only qualitatively analyzing. The reason of this lack of 
possibility of the determination of influence of the exhaust gas duct throttling on the exhaust gas 
flow. According to this the 3D model of the exhaust gas flow was prepared. The method applied in 
proposed model is Computational Fluid Dynamic (CFD). 


2. Laboratory research 


The chosen object of research is 3-cylinder, four-stroke, turbocharged, laboratory engine. The 
engine 1s loaded by a generator, electrically connected to the water resistance. During tests the 
engine was fuelled by diesel oil and operated at a constant speed, equal to 750rpm. The engine 
load and speed, parameters of the turbocharger, systems of cooling, fuelling, lubricating, and air 
exchange were measured. Pressure, temperature and humidity of air were recorded by laboratory 
equipment also. All mentioned results were recorded with a sampling time of 1 second. The 
scheme of the laboratory stand is presented in [4] and the engine parameters are presented in 
Tab.1. 

Tab. 1. The laboratory engine parameters 


Parameter Value | Unit 
Max. electric power 250| kW 
Rotational speed 750 | rpm 
Cylinder number 3 = 
Cylinder diameter 250 | mm 
Stroke 300 | mm 
Compression ratio 12,7 = 


The experimental study consists of 3 stages of 3 observations with simulations of different 
malfunctions of the exhaust gas duct. During each start of the observation, the engine was loaded 
to maximum load equal 250kW, and, after stabilizing the temperature of the exhaust gas behind 
the turbine, the engine operating parameters were recorded for 3 to 5 minutes. After this, the load 
of the engine was decreased by 10kW and, after stabilizing the temperature of the exhaust gas 
behind the turbine, the engine operating parameters were recorded again. The observation was 
continued with loads up to 50kW. The engine did not work with a load of 190kW due to resonance 
vibrations. 

Stages of experiment were set as follows: 

— the first stage during the operation of the engine assumed as “working properly", 
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— the second stage during the operation of the engine with the throttling of the cross section 
area of the exhaust gas duct by changing the barrier angle mounted in the exhaust gas duct 
behind the turbine by 21 degrees, 

— the third stage during the operation of the engine with the throttling of the cross section area 
of the exhaust gas duct by changing the barrier angle mounted in the exhaust gas duct 
behind the turbine by 71 degrees. 

The exhaust gas duct throttling simulation consisted of changing the angle of the throttling 

barrier mounted in the exhaust gas duct behind the turbine. Scheme of the exhaust gas duct and the 
prepared 3D model are presented in Fig.1. 


Exhaust gas 
—— 


Fig.1. Scheme of the exhaust gas duct and the 3D model 


3. CFD model 


Obtained laboratory research results and the geometric dimensions of the exhaust gas duct 
allow preparing the exhaust gas flow model behind the turbine. Throttling of the exhaust gas duct 
by the rotating plate causes the decrease of the cross section area of the exhaust gas duct and the 
changes in the flow characteristic, resulting the changes of the duct geometry. The result of this is 
intensity of turbulence phenomena in the exhaust gas flow. 

According to the mentioned purpose the modelling of the exhaust gas flow through the exhaust 
gas duct was conducted. Chosen method of modelling was Computational Fluid Dynamic (CFD) 
method. 

The used software to CFD conduction was AVL Fire. The exhaust gas duct construction was 
transferred and meshed to the pre-processor of the AVL Fire package. The cubic mesh with the 
minimum dimension equal 2,5 mm was chosen. The area of the rotating plate was re-meshed to the 
minimum dimension equal 1mm. The large dimensions of the exhaust gas duct and its axis 
symmetric construction, only the axis-symmetric half of the exhaust gas duct was modelled. The 
momentum and the continuity equations with k-zet-f model of turbulence [10], [14] were 
performed. The energy balance equations, the heat transfer phenomena and the gravity forces were 
neglected. 

The boundary conditions were selected as follows: 

— the exhaust gas duct inlet with the ambient temperature and the ambient pressure of the 
exhaust gas behind turbine, 

— the exhaust gas duct outlet with the ambient temperature of the exhaust gas behind the 
turbine and the environment pressure, 

— the axis-symmetry in the axis-symmetry surface. 

The maximum number of iterations was set on 3000 and the convergence criteria of the 
normalized residuals were set. The calculated result, as an extrapolate values of the cross section 
area of the exhaust gas duct outlet, is mass flow of the exhaust gas. 

Obtain the correct results of the calculation depends on specify the quantity of the Prandtl 
number and the density of the exhaust gas [6], [9]. The exhaust gas is mixture of the gaseous 
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combustion products, the water steam and not combusted in the engine cylinders air. The results of 
the exhaust gas analyze, obtained from [5] allows to specify the exhaust gas composition. 

Basis on the NIST Chemistry Book [11], about the parameters of the gaseous compounds of 
the exhaust gas in the obtained temperature and pressure ranges, the calculations of the Prandtl 


number and the density of the exhaust gas were conducted. 
250kKWa:IT_3000:Flow:Velocity[m/s] 


250kwa:IT_3000:Flow:velocity[m/s] 
148.59 


250KWIT_3000:Flow:velocity[m/s] 


Fig.2. Velocity in the exhaust gas duct for 250kW the engine load 


Mentioned calculations were conducted according to the weighted average by the following 
formula: 


i=l 
ре У ВЕ, (1) 


were: 
Pr; — the Prandtl number for the i-th gaseous compound, 


U, — the fraction of the i-th gaseous compound in the exhaust gas. 


The density of the exhaust gas was calculated accordingly to the (1) formula. 
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4. Results and discussion 


Mentioned actions allow to obtain the mass flow of the exhaust gas through the exhaust gas 
duct of the laboratory engine. The calculations for all considered loads of the engine and all angle 
positions of the throttling plate presented in [5] are done. 

On the Fig.2 the example of calculation results in the form of exhaust gas flow velocity for the 
engine load equal 250kW and different positions of the throttling plate are presented. Mentioned 
results are qualitatively representative for all considered loads of the engine. According to 
presented results horizontal position of the throttling plate disturbs the exhaust gas flow only a 
small degree. The flow is laminar. Only flow near the duct walls causes the changes on the flow 
direction. The 21? angle rotation of the throttling plate causes disturbance of the exhaust gas flow. 
It should be noted that the maximum velocity of the exhaust gas flow in the central point of the 
cross section decreases by 10%. Further increase of the exhaust gas duct throttling decreases the 
mass flow of the exhaust gas. The maximum velocity in this condition decreases by 50%. 
Presented results are qualitatively similar to presented in [1], [8], [13]. 

Fig.3 presents the results of calculations of the exhaust gas mass flow and measured parameters 
of the exhaust gas behind the turbine. According to presented results increase of the exhaust gas 
duct throttling causes the decrease of the mass flow of the exhaust gas. 


а 16 Exhaust gas flow TM Pressure drop РЕ 
14 Ы 
п u 1900 - " 
нги u 
— 12 u = 1700 4 
$ a u Е е * 
= 1 А = 1500 + 
5 Фо e. a M 
2 08 е ° . * 9 1300 
s . m properly а АС 
r e21 $ 1100 * M 
2 06 a + e? n 
h +71 900 | + MD AL. 
s 700- @ "E ga + у e a" m properly 
° 
NE” + KOŁACH I” ии ы 
02 à + * 500 | + nu e21 
+++ ы ++* * +71 
0 T T ++ т 
40 90 140 190 240 
Power [KW] 
с Exhaust gas temperature d * Specific fuel consumption 
aot t : * 550 | © 
+ + e? ag 
400 é PA + 
229 a" 500 Ы 
o + o z a + 
= te = ° 
© + a 2450 + 
2350 сов $ u" et 
5 И : H Г] m properly © 40) = .. m properly 
2 
5 +." °21 5 u е e21 
a " м, 2 m +71 350 u e+ : +71 
300 en "- * oe 
, ou ". 9% 55 + 
= 300 - "- фе 5 
u LEE ń 
250 T T 250 — T T т 
40 90 140 190 240 40 90 140 190 240 
Power [kW] Power [kW] 


Fig.3. Measured and calculated parameters of the exhaust gas duct 


Throttling of the exhaust gas duct decreases the pressure of the exhaust gas behind the turbine 
and increases its temperature behind the cylinders. Fig.3b and Fig.3c shows changes of the exhaust 
gas pressure and temperature behind the turbine. According to the presented results, the 
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temperature of the exhaust gas behind turbocharger increases with the increase of both exhaust 
duct throttling and the load of the engine. The exhaust duct throttling causes abnormalities in the 
combustion process. Abnormalities of combustion process cause a growth of the fuel consumption. 
Fig.3d shows specific fuel consumption measured for all considered states of the engine. 
Throttling of the exhaust gas duct by changing the barrier angle by 21deg. causes average increase 
of the fuel consumption by 3,796, and 8,2% for changing the barrier angle by 71deg respectively. 

The decrease of the differences of the exhaust gas pressure in the exhaust gas duct causes 
decrease of the turbocharger performance. The quantity of air delivered to the engine decreases. 
Simultaneously with this phenomena, the quantity of delivered fuel increases. The effect of this is 
combustion of the richer mixtures in the cylinders. 
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Fig.4. Approximation coefficients for exhaust gas duct flow 


The second goal of this research was approximation of the mass flow through the exhaust gas 
duct in accordance to different position of the throttling plate. Mentioned approximation may be 
useful tool to modelling the combustion processes in the engine cylinders in the different than 
laboratory conditions. 

Fig.3 presents the linear dependence between the exhaust gas mass flow and the angular 
position of the throttling plate of the engine in the constant rotational speed. According to this the 
approximation of this dependence is the linear function as follows: 


q=A-An+B, (2) 


where: 


q — the mass flow in [kg/s], 
An - the angle position of the throttling plate [?], 
A, B — coefficients of the linear function. 


Change of the engine output power causes changing the values of the mentioned coefficients 
“А” and “В”. Fig.4 presents the changes of the coefficients values in function of the engine power. 
The black dots presented the measured points. According to presented results the values of both 
“А” and “В” coefficients may to be approximated by a second order polynomial, according to the 
following formulas: 


A =а1.Р? +а2.Р a3, (3) 
B=b1-P?+b2-P+b3, (4) 
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Where: 


P — the power output of the laboratory engine [kW], 
al + a3, bl + b3 — the coefficients of the second order polynomials [-]. 


Presented analyse allows to select the values of values of coefficients from “al” to “ЪЗ” for the 
considered laboratory stand. Mentioned values are presented in Table 2. The results of calculation 
the “А” and “В” coefficients with use of values from the Tab.2 are presented in Fig.4 by 
continuous line. 


Tab.2. The coefficients values of the approximation function 


Coefficient al a2 a3 bl b2 b3 
Value [-] -3E-07 | 8E-05 | 0,0182 | 2E-05 | 0,0062 | 1,5124 


Presented results allow calculating the quantity of the exhaust gas mass flow for different 
values of the exhaust gas duct throttling in the range from the 0 to the 71? of the angular position 
of the throttling plate area without direct measurement. The maximum error of obtained 
approximation compared to measured values not excided 11%. 


5. Conclusions 


The main purpose of the presented manuscript is determination of the flow characteristic of the 
exhaust gas duct of the marine 4-stroke diesel engine. The simulation of the exhaust gas duct 
throttling by rotational throttling plate was done. The 3-dimensional flow model witch use of CFD 
method is applied. Obtained calculation results allow to determination of the exhaust gas mass 
flow for the simulated flow characteristics. The model of turbulence flow was taken into account. 
The gravity forces and the heat transfer phenomena were neglected. Obtained calculation results 
are qualitatively consistent with results obtained from literature. The analyze of the velocity 
distribution in the exhaust gas duct allows to conclusion that the changes of the angular position of 
the throttling plate causes significant disturbances in the exhaust gas flow. This causes the eddies 
creation and back flows. The result of this is the decrease of the exhaust gas flow. 

Additional purpose of the manuscript was approximation of the obtained results of the exhaust 
gas flow for different angular positions of the throttling plate. Obtained polynomial function may 
be useful tool to modeling the combustion process in the engine cylinders for the different flow 
characteristics of the exhaust gas duct. The calculation results allow to determination the mass 
flow of the exhaust gas with mean error equal 1196. The author's assessment is that obtained 
accuracy is not enough to use in the modeling practice. Therefore, further work is required in this 
area. May find the following a parametric assessment of the exhaust gas flow phenomena. 
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Abstract 


A portable engine test bench makes it possible to measure power and torque, acceleration of given speed, 
breaking distance, and distortions of the speedometer reading. The operation principle of а portable test bench is 
similar to that of a stationary load test bench with the only difference that the vehicle testing takes place in real road 
conditions, not in simulated ones, on the rollers of a stationary engine test bench. 

In this work, the operation principle of the portable engine test bench, the measurement method and exemplary 
results that have been obtained from carried out experimental test have been presented which provides the basis for 
further measurements of the vehicle operation real parameters in real not laboratory conditions. 


Keywords: measurement of engine power, portable engine test bench 
1. Introduction 


In their service instructions, car producers provide some engine parameters. This information is 
of great significance for the user, especially for the potential buyer of a given vehicle. From this 
data customers can draw fargoing conclusions. 

These parameters include: 

- engine capacity 

- torque 

- power 

As it is commonly known, the engine power means work of pressure of gases acting on the 
piston bottom, performed in a time unit. For multi cylinder engines it is the sum of the power of 
cylinders that make up the power unit. 

The driver is interested in useful power. This is power which can be transmitted to a receiver in 
any conditions of the engine operation [1, 3]. 

Rated power of an engine is guaranteed by the producer of the power unit for given conditions 
of operation. For passenger car engines the rated power equals maximal power, i.e. such that the 
engine can provide under constant load over a given time without the fear of exceeding 
mechanical allowable load or overheating. The maximal value of power is given in most of service 
instructions and technical data specifications for passenger cars. Both power and torque are 
variable in the function of the engine rotational speed. 
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A portable engine test bench makes it possible to measure power and torque, acceleration of 
any speed, breaking distance and distortions of the speedometer reading. The operation principle 
of the engine test bench is similar to a stationary load engine test bench with the difference that 
the measurement is performed in real conditions on the road not in conditions simulated on rollers. 

A portable engine test bench has many advantages [2]: 

- accuracy — with the assumption that the parameters are appropriate, measurement error is 
smaller than 1%. 

- the car is tested in its natural conditions, on the road — there is no problem with improper air 
supply, which makes the measurement accuracy higher than in a test house on a chassis 
test bench, 

- frequency of drawing data onto a chart: 60 times per a wheel revolution. Power and torque 
is provided with accuracy up to 0.1 KM and 0.1Nm. and the revolutions for which these 
parameters are obtained with accuracy of 1 revolution. . 

- program for control of the engine test bench enables manual entering of transmission 
values — it allows to eliminate errors which occur due to wrong calculations by programs 
for control of engine test bench transmissions. 

- program for control of the engine test bench has many interesting functions: ability ot 
compare up to 4 charts simultaneously, clear results in the form of tables and charts, full 
control of printouts etc. 

- correction of measurements according to DIN70020 ensures that neither the air temperature 
nor pressure affects fluctuations of power measurements. 

- there is a possibility to test a car with four wheel drive. 


Fig. 1. Image of portable engine test bench 
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2. Operation principle of portable engine test bench 


Like ша load engine test bench: the car, being in a given gear, is accelerated within an 
established range, then it is switched into the neutral and the rolling resistance is measured- with 
the difference that the sensor reading 60 impulses per one wheel revolution is fixed on the car, 
not at the rollers of a stationary test bench. 

Due to such an operation principle the device is characterized by high accuracy. For the 
purpose of high measurement accuracy it is necessary to be familiar with the vehicle weight, gear 
box ratio (it can be determined on the road), size of tires. It is also necessary to enter the air 
temperature and atmospheric pressure to make the given result recalculated with DIN norm 
obligatory in Europe. 

Since changes of atmospheric pressure, air humidity or its temperature affect the measured 
values — they are normalized and the obtained measurement results are converted into values 
corrected with normal parameters ( according to DIN 70020 — 25C?and 1000 HPa (does not 
account for humidity) and a special correction factor -power is estimated to change by 1% for 
temperature change by 5-6 degrees) engines tested according to American norm SAE do not have 
some accessories: air filter, alternator/ generator, fuel pump, clutch, radiator fan and exhaust 
system — therefore, the results according to SAE norm are higher by nearly 15-20 95. 


3. Description of the experiment 


The tests were performed with the use of Honda Accord VI from 2000 with engine 1.8 16V 
(136KM). Preparation for tests involved screwing a measuring disc and its revolution sensors to 
the vehicle wheel, connection of a cable for reading of the ignition system impulses (fig. 2)- 
measurement is performed in normal riding conditions and the results are transmitted to a portable 
computer which shows particular results of measurements.(fig. 3). 

In order to achieve proper results, such data as temperature, air pressure, size of tires, ratio of 
the power transmission system — were given. 


Fig.2. Manner of fixing a measurement disc with the sensor 
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Fig. 3. Image of Dynomet dk program window 


The measurements have proved that (fig.4) although the vehicle is 12 years old its power in 
real conditions does not differ from the power indicated by the producer. 
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Fig. 4. Exemplary chart from the tested vehicle power measurement 


During testing it was possible to observe the vehicle actual speed (fig.5) 
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Fig. 5. Chart of speed shows how the vehicle actual speed changes in te function of time. 
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Fig. 6. Measurement table 


In the above presented table 3 the first columns are: Time = (s) — Speed (km/h) and Distance Dist 
— distance covered (m). 


5Summary 


A device of this type can provide data on the car acceleration, its elasticity in particular gears 
within assigned values, power and torque in the function the engine revolutions and even the 
braking force. Accuracy of the power measurement on the wheels is increased by the value of air 
resistance — significant for high speeds, and it allows to establish what actual power is actually 
available and what is the maximal power that can be reached by the vehicle [4]. Undoubtedly, a 
disadvantage of such a device is approximate power of the engine (calculated by means of special 
algorithms which are designed to calculate approximate value of losses in the power transmission 
system). Moreover, when there is a sensor fixed on the wheel the measurement must be performed 
precisely along a straight line so that the wheel will not travel a distance different than the one 
that has been assumed. 
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This study aimed to carry out a preliminary measurement verification test, the development of 
reliable methods for measuring and obtaining data for further verification. 

While further research should be carried out verification of the above mentioned studies and 
research carried out by stationary dynamometer. Only the results of these studies will help to 
understand the real differences arising from the use of the vehicle in real time. 
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Abstract 


Problems connected with the circumferential distribution of compression ring pressure against the deformed 
cylinder wall during engine run were presented in this paper. In particular, division of the resultant force pressing the 
ring against the wall into component connected with ring own elasticity and the one originated from gas pressure was 
carried out and the influence of each force at individual phases of engine run was evaluated as well. The relationship 
allowing calculation of the circumferential ring pressure was also established. The paper contains figures and drafts 
obtained in a course of exemplary calculations of compression ring pressure concerning a running bulldozer engine. 
These figures show the areas of cylinder wall where so called light slots can appear. 

Conclusions resulting from the carried out simulation tests are discussed in Summary. The need for further 
investigations mentioned in presented study was shown as well. 


Keywords: combustion engine, piston ring, ring elastic pressure 


1. Introduction 


As soon as new cylinder liner is being installed a considerable deformation can occur caused 
by an erroneous installation in cylinder block. These deformations can grow to considerable size 
(of a few dozen micrometers [2]) during engine operation when the liner is being exposed to the 
influence of changing loads. Examples of such effect are presented in Fig. 1. 

The compression ring is being pressed to the cylinder wall with the radial forces that result 
from material own elasticity and gas pressure acting upon ring surfaces (see Fig. 2). While the 
elasticity force F, remains roughly constant during engine run, the other force F, changes itself. 
Investigations on the compression ring position in ring groove allow to conclude [1, 4] that the 
ring is pressed to the groove bottom side on the prevailing section of stroke which means that the 
acting gas pressure is almost the same as the pressure in combustion chamber (pa). Only on short 
sections of ring path the ring loses its contact with the bottom side of groove which causes that the 
F, force value depends above all on pressure рь in space below the ring. 

On the other hand the ring is separated from the cylinder wall with the Ру force resulting from 
the hydrodynamic pressure of oil film. The value of mean pressure is given by the equation [5] 
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qub, пу? 
Ps = Pat (Py — Р.И’, + h? W, + I? ns (1) 
where 
Pa pb — external pressures acting on the ring, 
u, v — ring velocities, axial and radial respectively, 
7] — viscosity of lubricating oil, 
hy — oil film minimum thickness, 
by — height of ring face covered with oil film, 
W; — load capacity coefficients, that allow to consider the effect of external pressures (p 
index), slide (u) and squeeze (v); their values depend on the profile of ring face. 
—PI 
—PII 
Se LA AAS Га үү 
г № 360% 07 gm м 
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Fig. 1. Deformations of cylinder sinusoidally 
loaded along its circumference; P I and P II — 


d 


Fig. 2. Sketch of compression ring moving 

relatively to cylinder wall; рь рь = external 

pressure, F,—force of ring own elasticity, F,— 
load of upper and lower part of rig ошл of oll layer, T friction beween ring 
respectively, 5 I and S II — form of bore before and oil layer. AA ring spee d, by—ring axial 
and after loading [3] height covered with oil film | 

If the value of oil film pressure was not enough to balance the forces pressing the ring against 
the wall the film rupture would occur and the ring would collaborate with wall in conditions of 
mixed or even boundary lubrication. In case of incomplete contact of ring face and cylinder wall 
the circumferential slots can occur which lead to the gas blow-by. This can cause the drop in 
engine power and an excessive oil consumption (which in turn leads to the increase in exhaust 
emission). Another task performed by rings, namely the heat conduction from piston to cylinder 
deteriorates as well. 

In his earlier papers author evaluated the ring pressure circumferential distribution against the 
deformed cylinder wall caused only by the ring own elasticity [7]. This pressure is additionally 
strengthened by the exhaust on the running engine which can lead to a different collaboration of 
these parts. Due to that the determination of combined effect of elasticity and gas force on pressure 
distribution and eventual definition of areas of ring and deformed wall contact was chosen as the 
goal of described tests. 

In order to realize the given aim the model tests were performed in which an extremely 
unfavorable case was selected, i.e. the lack of continuous oil film over the cylinder wall to 
eliminate the disturbance factors. Though this assumption is merely of theoretical character 
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(engine operation in such conditions would lead to a quick failure of collaborating parts) the 
obtained results allow to estimate value of some parameters indeterminable with other methods. 
One can mention here the forces loading radially the ring indispensable for full contact of ring and 
arbitrarily deformed cylinder wall. 


2. Ring to wall pressure from elastic and gas forces 


It was proved in [8] that the ring circumferential pressure distribution p,,(@) over the deformed 
cylinder wall can be determined according to the formula: 


Рь(Ф) = c [к Tu (z, | zlo) F2 zr) 27 (o). (2) 
pln 
where: 
E — Young's modulus, 
I — inertia moment of ring cross section, 
rm — radius of neutral layer of ring in cylinder, 
h, — ring axial height, 
z, — constant deformation of cylinder, 
2ь(ф) — deformation of cylinder wall, variable along the cylinder circumference (and its 
second and fourth derivative), 
К — characteristic coefficient of ring given with the formula 


Pm ity Tn 
a | ©) 


a 


Fig. 3. Sketch of a ring with marked 
dimensions and position of tangential force 
Е, апа resultant gas force F,- 


The value of K coefficient calculated according to formula (3) is just a reflection of ring 
dimensions and its material properties. However, taking into account the aim of investigation a 
trial was undertaken to such define the K coefficient that it would describe not only its elasticity 
but also the effect of exhaust gases. Using the following formula 


= т. 4 
ux (4) 


K 


that is another form of the formula (3) from [6], following relation has been obtained that allow to 
calculate a modified form of the К coefficient (because of its relation to crankshaft angle a further 
denoted as K-(a)): 
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K (2)- («Ee юка) ЕАР (а) т (5) 
E-I E-I 
where: 
Е, — tangential force fixed in the plane of neutral axis at ring free gap, 
F,-(a) — force fixed in the plane of neutral axis at ring free gap, equal to the gas force loading 
the ring, 
ре(о) — exhaust gas pressure loading the ring, 
a — crankshaft angle. 


Gas pressure loading the internal and working ring surface depends on the instantaneous 
pressure over (ра) and below (рь) the ring and on the position of ring in groove. For an initially 
adopted assumption about the lack of continuous oil film over the liner surface a linear drop in 
pressure along the ring face axial height was assumed. Due to that one can distinguish two cases 
for which the approximate value of gas pressure loading the ring is: 


Pela) =0.5-(p,(@) - p,(a)), (62) 
when the ring leans on the bottom side of groove and 
p,(2) 20.5 (p, (a) - p,(a)), (6b) 
when the ring leans on the upper side of groove. 
When the pressure over and below the ring is equal, the resultant pressure is zero, which means 
that both characteristic parameters аге the same, i.e. К = K,(a). 
3. Calculations of compression ring pressure on the deformed cylinder wall 
The presented below example concerns a trial of computational definition of compression ring 


wall pressure distribution against the deformed cylinder of earth moving machine engine (DT 466 
International Harvester). Technical data of this ring have been presented in Table 1. 


Tab. 1. Technical data of exemplary IC engine compression rings 


Ring 
Quantity 

cylinder diameter d [m] 0.109 
ring neutral radius rg [m] 0.0522 
axial height h, [m] 0.003 
radial thickness gp [m] 0.0046 
Young modulus E [Pa] 11210? 
mean pressure po [MPa] 0.178 
tangential force F, [N] 27.6 
stiffness EI [Nm] 2.65 
parameter K H 0.0286 


According to the primary goal of this study the further presented results concern an 
evaluation of the effect of gas forces on ring contact with a deformed cylinder wall. However, 
these deformations were limited to the one case, namely a circumferential line like in Fig. 1 
described with the formula 

Z, = A,cos(4p + д), (7) 
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where 44 and 6, are amplitude and phase shift of the Fourier series fourth harmonic, respectively 
(see Fig. 4). 
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Fig. 4. Course of cylinder circumferential Fig. 5. Course of K.(a) coefficient value vs. crank angle 
line for A,— 50 um and à; = x for full (1) and partial (2) engine load 


According to the analysis presented in Chapter 1 it was assumed that the ring is pressed against 
the bottom side of groove during the stroke of compression and expansion. The ring position 
during other strokes has been regarded as less important for the analyzed example because of far 
lower gas pressure (the radial position of ring depends above all on the ring own elasticity). 

Knowledge of engine constructional data as well as its operational conditions allowed to 
evaluate on the course of theoretical computations the pressure over and beneath the ring and 
eventually calculate the value of modified coefficient. corresponding to the consecutive positions 
of piston (only a fraction of these variations were presented in Fig. 5). The calculations of gas 
forces were carried out for engine full and partial load. 

Next, using Eqs. (2) and (7) the dependence (8) was defined that allows the calculation of ring 
wall pressure at an arbitrary point of cylinder circumference (9) and for a chosen crankshaft angle 
(a). 

Е-І 
Pn() = 7— 54K. (2) n, — 225. Ay cos4p + б,)]. (8) 


р Un 


Exemplary results presented further were achieved for the fourth harmonic amplitude which 
equals 10, 20, 50 and 100 um, respectively. 

The lines seen in Fig. 6 and 8 show the places on cylinder wall where the ring pressure falls to 
zero for a given engine load and amplitude of cylinder wall deformation (because of ring 
symmetry only its section from 0 to л/2 was shown). It was assumed that pressure equal or lower 
than zero means no contact of ring and wall. 

As it results from courses presented in Fig. 6 the area of ring face and cylinder wall contact 
increases along with the increase in ring wall pressure (caused by the increase of gases in 
combustion chamber). For example, during compression stroke for the cylinder deformation 
amplitude equal to 10 um the elastic ring pressure alone does not secure such contact. As soon as 
the gas pressure reaches the required magnitude which relates to the crank angle 260 degrees the 
contact is full. For the amplitude A, equal to 20 um the same effect is visible at 300 degrees CA 
while for 100 um the full circumferential contact appears close to TDC. During expansion stroke 
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the areas of full and partial contact have different shape because of different angular distribution of 
gas pressure. Reference [7] provides dependences that allow to calculate the amplitude boundary 
value of any harmonic which secures a positive ring wall pressure. On that basis the courses were 
drawn which allow to define the value of amplitudes for the analyzed example. 
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540 n op 
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Fig. 6. Course of boundary lines of positive ring wall pressure for ring moving during stroke of compression (a) and 
expansion (b) for engine full load and selected values of cylinder deformation; A4: 1 — 10 um, 2 — 20 um, 3 — 50 um 
and 4 — 100 um 
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Fig. 7. Fourth harmonic amplitude value 
20 vs. crank angle for engine full load at: 
1 — compression stroke, 2 —expansion 
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Fig. 8. Course of boundary lines of positive ring wall pressure for ring moving during stroke of compression (a) and 
expansion (b) for an engine partial load and selected values of cylinder deformation; A4: 1 10 um, 2 — 20 um, 3 — 50 
um and 4 — 100 um 
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Along with the decrease in engine load and reduction of gas forces the ring wall pressure also 
decreases which one can observe in the course of ring boundary pressure (Fig. 8). For example, for 
the amplitude of 100 ит and engine full load the ring touches the wall with its entire 
circumference close to the TDC. For the same amplitude of cylinder deformation and engine 
reduced load the light slot within a considerable angle of cylinder circumference emerges at the 
same area. 


4. Summary and conclusions 


Taking into consideration the results of theoretical studies and results of exemplary 
calculations following conclusions can be drawn: 

—  thering wall pressure caused by gas forces exceeds several times the pressure from ring own 
elasticity at strokes of compression and expansion, 

— the ring elastic pressure ensures the full wall contact only for minor deformations of cylinder 
wall (less than 10 um for the analyzed example), 

— for major deformation of cylinder wall ring pressure depends mainly on gas forces during 
strokes of compression and expansion close to the TDC in particular. 

It is worth mentioning that the presented studies concern the case of no lubrication. When the 
presence of lubricating oil is taken into account, the minor slots are filled with oil and gas blow-by 
would not happen. 

Another problem that should be the subject of eventual tests but was completely ignored in 
presented study is the analysis of stress in ring material during ring collaboration with deformed 
(worn) cylinder wall. Such tests should define the border values of cylinder deformations which 
would not lead to the ring failure. 
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Abstract 


The paper includes the results of the experimental study of the disjoining pressure in the thin oil films with 
anisotropic properties on the grey cast iron piston rings. It is shown that the disjoining pressure in the thin film, in a 
mode of self-regulation, can automatically balance the normal load between a pair “ring — cylinder liner”. The 
received data is used to increase the reliability of marine diesel engines by preventing from sudden failures as the 
result of piston ring breakage. 


Keywords: isotherm, pressure, oil film, anisotropy, ellipsometry, cylinder-piston group, diesel 


1. Introduction 


The organization of modern low speed diesel engines operating requires to improve the 
control of the friction processes of the mating parts of the cylinder-piston group (CPG) and to 
protect from the emergency situations. The breakage of forced piston rings of low speed diesel 
engines is the most common cause of functional failure. However, the cause of this phenomenon 
have not been studied properly yet. 

According to the literary analysis modern ships monitor the condition and functions of the 
individual essential CPG parts, including piston rings [1,2,11]. But these systems are used in the 
case of high speed of piston movement at hydrodynamic lubrication with isotropic properties. 
The piston rings breakage, in particular, can occur during passing through scavenge ports as a 
result of lubrication breakdown [11]. Thus, there is a need to develop the methods of technical 
condition identification at low speeds, while the ring passes through scavenge ports, at that the 
lubrication regime is not hydrodynamic, and the film has anisotropic properties. 

It is known that while the piston ring moves through the scavenge port of the cylinder liner 
the lubricant film drastically reduces, but the friction coefficient has a value u = 0,12 + 0,18, that 
indicates the absence of dry friction [11]. There is the disjoining pressure in a thin boundary 
layer of cylinder oil, which increases rapidly as the film thickness reduces at low speeds of the 
piston ring, due to the structuring of lubricant molecules [4,7,8,9]. 

The theoretical research of a ring movement process through the scavenge ports, in a case 
when a film with anisotropic properties separates parts, is not possible without the experimental 
value of the disjoining pressure in this film [9]. 

Thus, the lack of the properties research of thin lubricant films with the anisotropic properties 
on cast iron surfaces, of which the ring as well as cylinder liner are manufactured of, is a 
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deterrent for the improvement of low speed diesel engine reliability by preventing from sudden 
rings breakage, therefore such studies are relevant. 

The work objective is an experimental study of the disjoining pressure in thin oil films on the 
surface of the cast iron piston rings of low speed diesel vessels. 


2. Study of the disjoining pressure isotherm in the cylinder oil films 


This study was conducted using typical for the low speed diesel engine vessels cylinder oil 
ENERGOL CLO 50M, and the piston ring taken from the main engine MAN B & W 7S46MC- 
C. 

The metallographic analysis showed that the ring consists of grey cast iron; the base has a 
ferritic-pearlitic structure. The hardness of the ring is HB2770. The surface was polished to a 
value R,0,4+0,6, which corresponds to the surface condition of the piston ring and edges of the 
cylinder liner in service. 

The isotherms of disjoining pressure in a thin lubricant oil film on some steel surfaces, used 
in diesel engine manufacturing, was studied in the works [6.7]. The beginning of the film 
formation occurs at different thicknesses and pressures in the film. The film thickness and 
maximum disjoining pressure level depends on the chemical composition of the metal. 

The disjoining pressure isotherms study was conducted on an experimental ellipsometrical 
plant [5], such scheme is presented below. 


ly Sy 


M. эм - ellipsometric microscope 


Fig 1. The installation scheme of the disjoining pressure isotherms study in thin oil films on metal surfaces. 
1- cast iron sample; 2- steel cone;3- oil film; 4-groove;5- pressure vessel;6- hose 


The cell of the film was performed as follows: a steel cone with a hole of 1 mm diameter (2) 
was installed on the polished cast iron sample (1). The oil film (3) was formed in it and 
connected with the hose (6) to a pressure vessel (5), wherein there was the investigated cylinder 
oil. 

On a Gibbs' theory there exist transition layers at boundary of any adjacent phase, whose 
properties differ from those of the bulk phase. In the case of interfacial layers overlapping the 
hydrostatic pressure in the thin layer is different from the pressure of the bulk phase the part of 
which is the film, ie, Pascal's law is not implemented in the thin film[10,11] . The additional 
pressure providing the thermodynamic equilibrium of the film was called the disjoining pressure. 
It can be both positive and negative. The dependence of II — the isotherm of disjoining pressure 
is the thermodynamic characteristic of the thin film of liquid. 

To measure the disjoining pressure it is necessary to provide the mechanical equilibrium of 
the wetting film, by means of the external pressure. If the system is in thermodynamic 
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equilibrium, and the disjoining pressure is positive, then its value will be the low (negative) 
pressure produced in the conjugate bulk phase, by lowering the pressure vessel. Conversely, if 
the disjoining pressure is negative, its value will be excessive hydrostatic pressure. In both cases, 
the disjoining pressure Ia) will be equal to the difference between the pressure P, on the film 
surface and the pressure Po in the bulk phase [4,5,10]. 


P = pg (Ho. Hy), (1) 
where: 
p -oil density, 
Но - initial altitude of vessel lowering, 
H, - final altitude of vessel lowering, 
g- acceleration of gravity. 


Ho and H; values were determined using a micrometer device respectively, with accuracy 
AH == 0,1 mm, that caused the disjoining pressure error AP = + I Па. 

The experimental time to determine the equilibrium thickness for the wetting oil films with 
an opening of 1 mm diameter was 30 + 40 minutes. 

The position of the zero level Ho in the pressure vessel, corresponding to the time of film 
formation from the bulk phase, was determined according to the interference lines condition. It 
was set for the point at which the movement of convergence and divergence of the interference 
pattern was stopped. 

The film thickness В was measured using the ellipsometric microscope — М. The methods 
of the ellipsometric thickness measurement of wetting nonpolar oil films on conductive metal 
surfaces was developed by the authors and described in detail in the work [7]. The ellipse of the 
reflected light polarization is described by the ellipsometric angles ¥ and A, where tg Ч - is 
equal to the relative amplitudes change p- и S - a component, and A - the relative phase 


difference between them. The azimuth of analyzerA, , A, and polarizer Р, ‚Р, of light 
extinction were determined during the experiment. 

The thickness of layer was calculated using Drude equation [4], which connected the 
experimental parameters Ф and A and the optical characteristics of the reflected sample, defined 
with generalized Frenel coefficients Ry и Rs 


К 
iA P 
tgy et ==" 2) 
where: 
tg Y - relative amplitudes change, 
A- phase shift between the P and S components that arises during reflection, 
R,,R. - Fresnel generalized coefficients for the reflected light. 


The equitation (2) has the following form for the presented case of a homogeneous isotropic 
layer [4]: 


р’ 


-2ió -2ió 
tow: e^ — Rp, + Rye | Ее (3) 
EI 215 › 
1 F Ry, Re Rizs F Re 
2% А 
Here 9-7 hy n,dsin®, where d is the required layer thickness, and 
Ryu, Ria. Ras, Ras; are defined accordingly: 
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The task is complicated by the need to take into account the effects of the electromagnetic 
waves weakening in the metal substrate. In general case, the equation (3) is formulated for the 
inhomogeneous wave and can be solved by the introduction of complex refractive index, taking 
into account the damping effects. In the presented case 73 > N, = п; - ік; and n, k- are the 
coefficients of substrate refraction and absorption. The following values n = 3,9, К = 6.96. are for 
cast iron. 

Thus, the equation (3) with (4), is a complex expression, and hence, the simultaneous 
equations is required to be solved to find the layer thickness: 


(R, + Ке ?®)(1 +Rp.R,e *") 
‚е^ р p 
Re| tgy-e (1 + ЕК (Ra. c Rue") И 
Iml tey. e” (Ra + Ке? | 4 BS) D 
SV (1 + Boa Вы E Rye?) 


The experimental isotherm of disjoining pressure of the oil film on the cast iron surface with 
roughness (R,0,4), obtained at a temperature of 295 ? K is shown in Fig. 2. The isotherm is close 
to an exponential form that corresponds to the works 


[6]. 
A 
P = Pe Й (6) 


А - Hamaker constant 
The isotherms in the oil films on cast iron corresponds to the area of positive values of the 
disjoining pressure IIs > 0. They are of the decreasing nature. From this it follows that the 
disjoining pressure in the thin films, in the mode of self-regulation, can automatically balance the 


normal load, between the pair "ring — cylinder sleeve". Each new value of pressure has new 
equilibrium thickness of the film. 
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Fig. 2. The disjoining pressure change in a thin film of lubricant ENERGOL CIO 50M on the cast iron piston 
ring, depending on the film thickness 
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The maximum value of the disjoining pressure was determined from the graph constructed in 
the semi-logarithmic coordinates of dependence InP (h) on the film thickness h. The maximum 
value of the disjoining pressures for the cast iron is equal to Il =1400001I1a . 

The distinctive feature of the oil films with anisotropic properties on cast iron is the small 
amount of change in the film thickness from the beginning of its formation to the minimum at 
which the pressure in the film begins to increase sharply. For the films on the cast iron the 
thickness is 10 + 12 nm. In comparison with the thickness of the films on steels (Fig. 3) [6,7] we 
can see that on steel 35 XMA,of which the piston heads and MAN B & W were manufactured, it 
is equal to 126 nm, which is an order of magnitude greater than that of cast iron. 

The value of the disjoining pressure is the sum of the molecular component Птһ), at low 
film thicknesses and structural Ts): 


Пъ= IImgy* IIS). 


According to the obtained results the chemical composition of the substrate has a significant 
effect of on the change in the degree of orientation molecules order as the film thickness 
changes, that determines the amount of the structural constituent of the disjoining pressure Is). 
Obviously, that such a difference in the obtained values of II) can be explained by the different 
contributions of the molecular component of the disjoining pressure, which is essential for small 
film thicknesses. Previously, the effect of the molecular components on the disjoining pressure 
value was defined for polar liquid film on dielectrics [3,5]. 


Ila 30 


120 140 160 180 200 220 240 260 h*10?w 


Fig.3. The dependence of the disjoining pressure in the oil film ENEKOL GLO 50M on its thickness on steel 
SH15 


The analysis of the received values of disjoining pressure shows that the use of Is) as the 
exponent is not always correct, and the real dependence of structural constituent on the film 
thickness is more complicated. 

The analysis of the photomicrographs of the thin oil film on the surface of iron shows that the 
film is not flat by its nature. (Fig. 4) 

The distribution of the interference lines demonstrates the formation of the film with the 
anisotropic properties in the center of the cell and the increase of the film area as it is growing 
thinner (Fig 4.b, c, d, e). An interference pattern is formed by the graphite inclusions (Fig.4f). 

The control of lubricating film with anisotropic properties, determined by the disjoining 
pressure, make it possible to compensate the normal load on friction zone, therefore energy loss 
and the amount of conjugated surfaces deterioration can be reduced. 
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The obtained results enable to improve the identification methods of piston ring working 
capacity according to the lubricant state in order to prevent an emergency on board. 


e) 


Fig. 4.The photomicrographs of the oil film on cast iron: formation (a), the minimum film thickness (f) 


Conclusions 

1. For the first time the isotherm of the disjoining pressure in the cylinder oil films 
on cast iron was experimentally obtained. It corresponds to the positive pressure region 
Пу. There is force in a thin cylinder oil film directed to the opposite side of the ring 
pressure on the port of cylinder liner. 

2. Itis determined that the isotherm of disjoining pressure in the cylinder oil films on 
iron and steel surfaces are of the decreasing nature. It ensures the process of self- 
regulation of piston rings pressure on the cylinder liner. Each new pressure value has a 
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new equilibrium thickness. The process of self-regulation is automatic and it does not 
require the human intervention. 

3. The comparison of the experimental isotherm of the disjoining pressure in the oil 
film on steel and iron surfaces have the following results: 

e the film formation on steel and cast iron occur at different film pressure and 
thickness. The film formation on steel occurs at interval of 240 + 280 nm, and with a 
much smaller thickness of 163 + 165 nm on grey cast iron; 

e the film thickness from the beginning of its formation and before the sharp 
pressure increase is 130 nm for steel. It is higher than on cast iron, that is only 10 + 12 
nm; 

e the minimum film thickness on the iron surface where the disjoining pressure has 
a maximum value is in the range of 155-160 nm; 

e the maximum disjoining pressure in the cylinder oil film was determined. It is 
equal to П;=0,14 МПа; 

The further development of the piston rings movement process study, in the presence of thin 
oil films with the anisotropic properties, using the obtained data of the disjoining pressure on 
cast iron, will enable to improve the reliability of marine low speed diesel engines. Identification 
methods of piston rings technical condition and prevention from the sudden failures in the result 
of the piston rings damage will be developed. 
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Abstract 


Low calorific gaseous fuels are becoming increasingly used as alternative fuels for feeding power equipment of 
all kinds, including combustion engines. The increased interest in these fuels results, above all, from support for the 
development of this type of power engineering by governmental and international institutions, aiming to increase the 
renewable energy share in the overall energy balance. Fuels of this type are obtained by processing different kinds of 
biomass and natural raw materials and can also be recovered from natural processes occurring in nature. Due to the 
diverse technologies for obtaining these fuels, their chemical composition is variable, which significantly affects the 
efficiency of their use as fuels for combustion engines. When engines are fed with gaseous fuels, combustion chamber 
filling conditions change considerably due to the much higher gaseous fuel volume compared to liquid fuels, which 
significantly affects the engine’s performance. 

This paper presents the effect of individual combustible components contained in different low calorific gaseous 
fuels on changes in the degree of engine combustion chamber air filling and on the calorific value of the produced 
combustible air-gaseous fuel mixture. 


Keywords: environmental protection, low calorific gaseous fuels, calorific value, excess air coefficient, combustible 
mixture 


Introduction 


Growing social awareness and concern for the environment and the increasingly stricter 
standards for reducing greenhouse gas and toxic compound atmospheric emission limits have 
prompted a search for new alternative energy sources [1, 6]. Another important factor causing an 
increased interest in renewable energy sources is the need to increase the share of energy produced 
from alternative energy sources as required by EU legislation. It is also worth stressing that the 
heightened interest in alternative fuels is favoured by shrinking resources of fossil fuels, 
particularly petroleum [4, 5, 6]. 

Low calorific gaseous fuels of different types are potential energy sources which can be 
obtained by processing different kinds of biological products and can also form spontaneously as a 
result of natural processes occurring, among others, in dumping sites, sewage treatment plants and 
animal farms [1, 4, 5]. 

The composition of such fuels is not constant and is closely dependent not only on the raw 
materials from which they are produced but also on the technology for their production and 
treatment. The combustible compounds making up fuels of this type include, above all, methane 


89 


(CH4), carbon monoxide (CO) and hydrogen (H2). The basic properties of the listed combustible 
components of gaseous fuels are presented in Table 1. 

The percentage of individual components in gaseous fuels is not constant and, as has been 
mentioned, is conditioned by many factors. The main low calorific gaseous fuels most often used 
currently include [1, 2, 5, 7, 8, 9]: 

- agricultural biogas — obtained by methane fermentation of plant raw materials or animal 
excrement; 

— biogas from sewage treatment plants — gas recovered in sewage treatment plants in sewage and 
water treatment; 

— landfill gases — forming by decomposition of organic substances contained in municipal waste 
after disposal; 

— fossil gases — gases obtained in hard coal mine demethanization; 

— synthesis gases (syngas) — obtained in gasification of solid fuels, petroleum refining residues, 
steam methane reforming or partial methane oxidation with oxygen. 

Table 2 presents the percentages of individual combustible components in the most common 
low calorific gaseous fuels. 


Table 1. Basic properties of combustible components of low calorific gaseous fuels 


Gas с | gay | Adema e 
[к/т] | Umm] [K] 
Hydrogen (H5) 10,78 0,0899 2,38 807 
Carbon monoxide (CO) 12,6 1,25 2,38 881 
Methane (CH4) 35,89 0,717 9,54 923 
Ethane (C2H6) 63,77 1,34 16,7 793 
Propane (СзНз) 91,28 1,97 23,8 783 
Butane (C4H;0) 122,6 2,33 31 692 


Table 2. Volumetric content of individual components in gaseous fuels [2, 4, 6, 8] 


Gaseous fuel e ii in к > R 
[^o] [^o] [7o] [^o] [7o] [7o] 
Agricultural biogas - - 45-75 | 25-55 0-5 0-3 
Biogas from sewage treatment - - 57-62 | 33-38 3-8 0-1 
plants 
Landfill gas 0-15 - 37-67 | 24-40 | 10-25 | 0-45 
Mine gases - - 40-80 | 8-15 | 20-40 | 0-10 
Synthesis gases 3-0-70 | 0-60 0-10 5-30 0-2 - 


Because of the quite large differences possible in the chemical composition of such gaseous 
fuels, the combustion conditions for such fuel in the engine change significantly, which causes 
significant changes in the obtainable engine power, the amount of heat carried away by the 
engine's cooling system and heat carried away by exhaust gas. 
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Characteristic parameters of the fuel-air mixture 


The power obtained from a combustion engine depends primarily on the calorific value and the 
amount of fuel burnt in the combustion chamber. The amount of burnable fuel is conditioned by 
the amount of oxygen available in the engine chamber, which is necessary for burning combustible 
components contained in the fuel. The vast majority of the oxygen supplied to the engine's 
combustion chamber comes from the air and only a small part of oxygen can be supplied with fuel 
(for oxygen-containing fuels). 

Therefore, the power obtainable by the engine is limited by the amount of air (oxygen) induced 
into the combustion chamber needed for complete fuel combustion. For liquid fuels, the theoretical 
amount of induced oxygen results from its content in the air and the maximum combustion 
chamber volume. In this case, the volume of supplied liquid fuel is minute and can be ignored 
practically without error. In reality, the induced fresh air volume depends on the engine filling 
efficiency, which is most affected by the structure of ducts feeding fresh load to the combustion 
chamber, the efficiency of cleaning exhaust gas residues from the previous work cycle of the 
combustion chamber, engine speed as well as the temperature in the combustion chamber during 
the suction stroke. In real combustion engines, cylinder filling efficiency ranges from 0.7—0.85 and 
decreases with rising engine speed, due to shortening suction stroke length. 

While for liquid fuels the volume of fuel fed to the combustion chamber can be ignored, for 
gaseous fuels, containing not only combustible compounds but also substantial percentages of 
non-combustible compounds such as carbon dioxide (CO2) and nitrogen (№), the volume of fuel 
fed to the engine chamber is considerable (except for solutions with gas injection directly into the 
combustion chamber). Therefore, this significantly reduces the volume of air induced into the 
combustion chamber for fuel combustion which, in turn, changes the engine performance. 

As mentioned before, low calorific gaseous fuels contain three basic combustible components: 
methane (CH4), hydrogen (Н2) and carbon monoxide (CO). Assuming that the whole amount of 
fed fuel is burnt, the occurring combustion reactions for these compounds can be written as: 


CH, + 20; = CO; + 2H;0 (1) 
2H; + О =2Н›О (2) 
2СО + О = 2CO; (3) 


Considering that 21% of oxygen is in the air by volume, it can be stated that the theoretical 
amount of air needed for the combustion of a fuel which is a mixture of these three gases is: 


L = оон, +С0)+2СН.] [n air / m? gaseous fuel] (4) 


For gaseous fuels containing oxygen, this amount of oxygen should be taken into account, 
therefore, the above formula assumes the form: 


1 А 
І, = он, + СО)+ 2CH , — O, | [m? air / m? gaseous fuel] (5) 
To determine the volume of air induced into the combustion chamber with the air-gas mixture, 
the degree of combustion chamber air filling у» can be defined. With the assumption that the 
cylinder filling efficiency is 7,= 1, the exhaust gas removal from the combustion chamber is 
complete and an air-gas mixture has formed outside the combustion chamber, this ratio can be 
written as: 


ү, VaV, 
Yair =. = " c (6) 
Vy, y, 
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where: 

Иш — air volume in the air-gas mixture; 

Ven — combustion chamber volume; 

Vjuel — gaseous fuel volume. 

Effect of gaseous fuel composition on the degree of combustion chamber air filling and the 
calorific value of the mixture 


The variable composition of gaseous fuels significantly affects both the degree of combustion 
chamber air filling and the calorific value of the produced combustible mixture. Fig. 1 presents 
changes in the degree of engine combustion chamber air filling for basic combustible components 
of gaseous fuels as a function of the excess air coefficient A. The presented changes clearly show 
that for feeding with gaseous fuels the volume of induced air is much lower than when the engine 
is fed with liquid fuels. For a stoichiometric mixture (A = 1), the volume of air induced into the 
combustion chamber is 1096 lower for feeding with methane and 3096 lower for feeding with 
hydrogen or carbon monoxide. Fig. 2 presents changes in the calorific value of the combustible 
mixture produced from individual gaseous fuel components depending on A. For comparison, this 
figure also contains changes in the calorific value of propane (C3Hs) and butane (C4Hi0), the main 
LPG components. 
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Fig. І. Degree of combustion chamber air filling for feeding with CH,, Н and CO as a function of A 
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Fig. 2. Change in the calorific value of the air-gas mixture depending оп A for different gaseous fuel components 


Fig. 3 presents changes in the degree of combustion chamber air filling Yair and changes in the 
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calorific value of the air-gaseous fuel mixture containing only biogas (e.g. methane) depending on 
the percentage of CH; in fuel with the excess air coefficient A = 1. It can be concluded based on 
changes in combustion chamber air filling that a rise in the percentage of methane in fuel 
significantly affects the use of the combustion chamber. For fuels with a ca. 40% methane content 
the degree of engine combustion chamber volume use is only 65%. 
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Fig. 3. Change in the degree of combustion chamber air filling and the calorific value of the air-gaseous fuel mixture 
for fuels with different methane percentages for A = 1 


Fig. 4 presents changes in the degree of combustion chamber air filling Yair with A= 1 for 
gaseous fuels containing two combustible components, e.g. landfill gases. Fig. 5 presents changes 
in the calorific value of the air-landfill gas combustible mixture for different concentrations of 


individual combustible components depending on the excess air coefficient A. 
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and Н; percentages for A = 1 
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— — CH4- 60%, Н2-5% 
— — CH4- 50%, Н2-10% 
— — CH4- 40%, H2-5% 
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Fig. 5. Change in the calorific value of the air-landfill gas mixture depending on A for different percentages of 
combustible components 


For the combustion of synthesis gases, obtained by gasification of biomaterials (e.g. biomass, 
wood), coal gasification or steam methane reforming, oxygen contained in this fuel should also be 
taken into account. In fuels of this type, the oxygen volume fraction can reach even 2096, which 
significantly changes fuel combustion conditions. Fig. 6 presents changes in the combustion 
chamber air filling ratio for selected synthesis gases [3] and Fig. 7 presents changes in the calorific 
value of the air-synthesis gas mixture depending on the pure oxygen content in this fuel for A= 1. 
The presented changes show that oxygen contained in fuel considerably increases the value of heat 
obtainable from the combustion of a specific portion of the produced air-gaseous fuel mixture. 
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Fig. 6. Effect of the oxygen content in the synthesis gas on the value of the combustion chamber air filling ratio 
maintaining A = 1 
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Fig. 7. Effect of the oxygen content in the synthesis gas on the calorific value of the combustible mixture for A = 1 
Conclusions 


In summarizing the results on the effect of the composition of low calorific gaseous fuels on 
combustion chamber air filling and the calorific value of the obtained air-gaseous fuel mixture 
with the combustible mixture formed outside the combustion chamber, it must be concluded that: 

- feeding the engine with gaseous fuels significantly reduces combustion chamber fresh air 
filling, which limits the burnable amount of fuel. When the engine 15 fed with pure methane 
with A = 1, the combustion chamber volume fillable by air decreases by around 10%. When the 
engine 15 fed with pure carbon monoxide or methane, the volume of induced air under the same 
conditions decreases by 3096. This reduces the amount of burnable fuel, which limits the 
amount of heat released in the combustion chamber; 

- when the engine is fed with low calorific fuels containing ballast compounds besides 
combustible components, the amount of heat obtainable in the combustion chamber is still 
lower, which causes a further decrease not only in the obtainable engine power, but also causes 
increased losses of heat absorbed by gases not participating in the combustion reaction; 

— for fuels containing oxygen, e.g. synthesis gases, the oxygen contained in fuel increases the 
calorific value of the formed combustible mixture along with an increasing oxygen percentage 
in the fuel; 

— it should, however, be stressed that despite the above limitations, low calorific gaseous fuels 
positively affect the share of renewable fuels in the overall energy balance. It is also important 
to stress that the use of these fuels for energy purposes reduces the emission of methane 
released spontaneously from, among others, dumping sites, which is ca. 25 times more harmful 
to the ozone layer than carbon dioxide [3]. 
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Abstract 


In the paper a physical model of certain processes taking place when marine diesel engine is driving 
a synchronous generator is presented. This physical model will be a basis for a mathematical model realization 
in a form of equations describing the relationships between the designated parameters of the diesel engine and 
the generator. The method of theoretical and operational investigations which enables identification of the 
technical condition of the set is presented. Systematic collecting the results in the database will allow to refine 
the model, and to conclude about the technical condition of real objects. 


Keywords: diagnostic, marine diesel generator 
NOMENCLATURE 


PARAMETERS 

- force 

- current 

- mass polar moment of inertia 
- enthalpy 

- torque 

- mass 

- mass flow rate 
- pressure 

- gas constant 

- heat flux 

- area 


Uic md HB deeem m 
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- temperature 

- voltage 

- linear velocity 

- volume 

- generalized losses 
- isentropic exponent 
- angular velocity 


82044 сч 


ABBREVIATIONS AND INDEXES 


cyl - cylinder 

8 - Bas 

K - connecting rod 

L - phase 

oK - connecting rod response 

oMP - engines machinery response 

oP - generator response 

oWK - crankshaft response 

p - back pressure 

P - generator 

pk - exhaust gases in the manifold 
sp - exhaust gases in the outlet channel 
t - piston 

wk - crankshaft 

* - total value of parameter 
1,2,3...n - subsequent cylinder, phase, etc. 


1. INTRODUCTION 


One of the maintain problems of the generating set consisting of the marine diesel engine 
and the synchronous generator is to determine their technical condition. The generators on the 
ship are frequently driven by a combustion engines that are not equipped with indicator 
valves. Such engines with limited monitoring susceptibility, are operated according to so- 
called overhaul life strategy. The current trend is to implement economically advantageous 
operating strategy according to the technical condition of the ship's mechanisms and 
appliances. Such attempts have also been made with respect to diesel-electric generating sets. 
Technical condition of selected elements of the structure design based on endoscopic 
examination [4] and vibration testing was evaluated. An alternative method currently 
developed by the authors of the paper is to determine the technical state of the ship generating 
set by means of testing pressure pulsations in the exhaust gas channels, phase-to-phase 
generator voltages, and vibrations of selected parts of the engine and the driven 
generator [2, 3]. 


2. PLAN AND METHODOLGY OF RESEARCH 


The first stage of research directed on developing methods for evaluation technical 
condition of the generating set of the ship power station was to develop a research plan. 
Similar plans were worked out in previous works by the authors, for example [2, 3, 5, 8] in 
the form of block diagram presenting successive stages of the investigation. The diagram 
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shows the interactions between particular stages of work. In general, the starting point for the 
development of the plan is a real object, in this case a diesel- electric generating set consisting 
of a combustion engine type SW 400 and synchronous generator GCPf-94c/1. Basic technical 
data of the investigated engine are shown in Table 1. 


Table 1. Technical characteristic of diesel engine type SW 400 [9] 


Rated output 54,06 kW | Injection type Direct 

Rated crankshaft 1500rpm | Injection order 1-5-3-6-2-4 

Speed 

Number of cylinders 6 Injection pressure 16,18 — 16,67 

MPa 

Piston stroke 120,65 mm | Specific fuel 190 g/kWh 
consumption 

Cylinder bore 107,19 mm | Intake valve opening 10° before the 
commencement TDC 


General view and schematic diagram of the generating set is shown in Figure 1. 


Fig.1. General view and schematic diagram of the diesel-electric generating set being the object of the research 
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In the proposed research plan the real object is a base for developing a physical model of 
the selected processes occurring in diesel-electric generating set. In this physical model the 
possibility to simulate damages to the combustion engine and synchronous generator by 
means of changing values of selected structural parameters was taken into account [6]. The 
physical model is a base for preparing a mathematical model in a form of set of equations 
describing the interdependence between selected mechanical, electrical, and gas dynamical 
parameters. 


The equations of the mathematical model recorded as a computer programme will allow 
for conducting the model studies of the effects of simulated damage on the waveforms of the 
generating set analyzed energetic parameters. At the same time we are going to conduct the 
research on the real object, generating set. There will be introduced the same known and 
recognizable damage to the model and the real object. The results obtained from simulation 
studies and research conducted on the real object will be compared in the time domain and 
frequency domain. Following the results of the comparative analysis the adequacy of the 
physical model and based on it the mathematical model will be evaluated. If the simulation 
results will differ to a greater extent than assumed by the results of a real object, the physical 
model will be modified to increase its relevance. The results of all studies, both the model and 
simulation will be collected in a database technical condition — symptom. Creation of such 
a comprehensive database will help in the future to identify the technical state of the 
generating set and the proper diagnosis. The proposed flow chart of the investigation is shown 
in Figure 2 [1]. 


Fig.2. Diagram of realization the diagnostics research 
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3. A PHYSICAL MODEL 


The proposed physical model is a certain conception of theoretical representation 
of phenomena and processes which occur in considered real object — generating set. This 
model includes a set of assumptions that define a simplified image of the investigated object, 
and takes into account whole physical quantitative and qualitative relationships. The degree 
of model simplification results from both objective and subjective factors. The objective 
factors include the state of the art — the possibility to identify all the phenomena and 
processes, and competency to describe them by means of appropriate relationships. 
The subjective factors are the needs and technical conditions for the model analysis of the 
given object (generating set). In general, only selected stages of considered processes are 
modeled, with far-reaching simplification [6]. 


In this model, the physical input parameters are: 


e gas pressure in engine cylinders ру = f (a), 

e cylinder volume Vey, = f (a), 

e gas parameters inside cylinders, such as gas constant Кс = f (a) and isentropic 
exponent Key, = f (a), 

e engine load expressed by generator voltage and current О, І = f (a). 


АП these parameters are a function of the engine crank angle a. The next group of input 
parameters are the design parameters of the structure of the generating set (diesel engine, 
synchronous generator, coupling). The main of them are the piston diameter, connecting rod 
length, and crank arm length. Another group of parameters are connecting rod, piston, piston 
pin masses distribution, rotating elements moments of inertia (engine crankshaft, clutch, rotor 
of the generator, driven mechanisms etc.). 


Output parameters are: 


e exhaust gas pressure in the outlet channel pp, = f (x, T), 
e exhaust gas temperature in the outlet channel Tj, = f (x, т), 
• speed of movement of peak amplitude exhaust gas pressure wave Wp, = f (x, t). 


All these exhaust gas parameters in the outlet channel are a function of time and 
the channel length (coordinate cross-section in which they are measured). The last group 
of the output parameters are those measured on the synchronous generator driven by the 
diesel engine. In the presented model instant values of phase-to-phase voltages as time 
functions are taken into account. The proposed physical model is shown in Figure 3. 


In this model the engine cylinder were treated as zero-dimensional objects 
ie. independent variable for them is only time of the process. Engine cylinders are able 
to store a mass of thermodynamic medium marked as m,y,. The input parameters relating 
to the engine cylinders are pressures inside cylinders, cylinder volumes and parameters 
of working medium inside cylinders such as individual gas constant and isentropic exponent. 
In the model of the engine cylinder section mass of the piston with connecting rod was taken 
into account, and marked as тт, and considered as carrying out reciprocating motion only. 


Engine cylinders output thermodynamic parameters being the input parameters for the 
exhaust channel are pressure рг, and temperature Тсу, of the exhaust gases, their mass flux 
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Mey, enthalpy flux їс, and speed of moving working medium ил. Additionally, the model 
takes into account influence of the gas pressure from the outlet channel on cylinder sections 
рр. Exhaust gas output channel is modeled as a one-dimensional object, which means that the 
independent variables are time and channel length determining the position of control cross- 
sections in which exhaust gas measurements are taken. The physical model takes into account 
dynamics of elements reciprocating and those in rotating motion. For this purpose weight 
distribution of selected structural elements of the engine and generator were taken into 
consideration. 


This refers mainly to the piston with piston pin and rings involved in reciprocating 

motion, connecting rod that moves in a mixed way — its one part is reciprocating and the other 
is rotating, and engine crankshaft with a flywheel that rotate together with the rotor of the 
generator. The complicated character of movement of the connecting rod can be examined as 
a sum of a motion of reciprocating mass and a motion of rotating mass, calculated by means 
of the “weighing system" [7]. In the model gas forces acting on the piston F,and connecting 
rod forces Рок were considered as well as moments of inertia of the piston-crank mechanism 
Jy. Engine crankshaft 15 another element of the model. Turning moment acting upon 
the crankshaft, originating from gas forces, is a product of these forces and the length 
of the crankshaft crank My. In a similar way as in the case of the piston-connecting rod set, 
in the assembly connecting rod-crankshaft the moment Moy; from the crankshaft acting 
on the connecting rod was taken into account . It is assumed in the model that all the masses 
and moments of inertia of mechanisms driven by the engine are treated globally. 
Torque derived from the crankshaft is transmitted to the generator rotor Mwg. The generator 
influence on the engine is the load torque Мор. It was assumed in the model of the generator 
that the only rotating mass is the mass of the rotor тр. The synchronous generator unit 
produces three-phase alternating voltage with a frequency of 50Hz which meets double 
frequency of rotating speed of the engine crankshaft. Electric load of the generator brings 
about passage of current I, at the assumption that phase-to-phase voltage is U,. As the 
generator load is resistive, power can be calculated as the product of voltage and current 
because power factor is equal to 1. Consumption of power produced by the synchronous 
generator is the source of the diesel engine load torque Мор. 


There are feedback loops in the model because increasing of the generator load results in 
decreasing of electric voltages in phases Lı, [2 and 1з what in turn increases generator 
excitation current. This process is controlled by AVG (Automatic Voltage Regulator) which 
senses the generator output voltage and acts to alter the field current to maintain the voltage at 
its set value Vi, V», and Уз. As a result of such a regulation the generator voltage increases 
or decreases to +2,5% (or better), so changing the angular speed of the diesel crankshaft о. 
The angular speed is an input signal to the speed governor of the diesel engine that, by means 
of changing the fuel charge, keeps the speed at the set value. This speed is an input to the 
speed governor of a piston engine, which is a drop in the value of increasing the flow of fuel 
supplied to the engine cylinders thus increasing the angular velocity of the crankshaft to the 
specified value. 
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Fig.3. Physical model of energetic processes taking place in the generating set diesel engine — synchronous 
generator 


4. SUMMARY 


The presented above research plan and physical model of the marine generating 
set should allow to work out a non-invasive method for evaluation the technical condition 
of a diesel electric set, and further a diagnostic methodology, when the engine is not 
equipped with indicator valves. It is estimated that diagnostic inference there will be possible 
about combustion process in engine cylinders. Technical state of the engine fuel system, its 
regulation and the tightness of the set piston-piston rings-cylinder liner, and the condition 
of the timing system, have the greatest influence on it. Moreover, it will be possible to assess 
the technical condition of the electrical elements of the synchronous generator and its 
bearings. 
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Abstract 


Transient states are an important part of the spectrum of engine loads, especially the traction engines. In the 
case of marine diesel engines, transient states are of particular importance in reducing the analysis of motion 
units for special areas and maneuvering in port, the participation of transient states in the load spectrum 
significantly increases, also, the emission of toxic compounds from this period increases proportionally. The 
factors which determine the value of the emission are the forces shaping transient states and the technical 
condition of the engine itself. This paper presents a description of transient states using multi-equation models, 
and the analysis of their relevance. It also presents a comparison of toxic compounds concentration at modified 
angles of fuel injection advance. 


Keywords: diagnostic, theory of experiments, marine diesel engine, exhaust gas toxicity, multi-equation models 
1. Introduction 


Transient states are exceptional marine diesel engine operating conditions. They arise in 
the absence of thermodynamic equilibrium in the engine cylinders and are an important part 
of the engine load spectrum, especially of traction engines, thereby without affecting the 
emission of toxic compounds. Engine research in this area is forced because of homologation, 
where the main problem comes down to the optimization of the combustion course with 
variable engine load described even through urban driving tests. 

In the case of marine diesel propulsion, the importance of transient states, in the above 
sense, is less prominent because of the relatively small proportion of transients in the engine 
load spectrum. If, however, such an analysis is subjected to the movement of individuals in 
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specific areas or maneuvering in port, the proportion of transients in the engine load spectrum 
grows significantly and is worthy of special consideration. Proportionately to this growth 
increases the emission of toxic compounds, caused by the impact of those states. This should 
be explained by the fact that transients interfere with cylinder thermodynamic equilibrium, 
which occurs during the fixed charges. This interrupts the combustion process by causing 
temporary changes, primarily to the stream of fresh charge of the cylinder, but also the 
amount of fuel delivered. Thus, the air-fuel ratio changes temporarily, which results in the 
changes in air excess ratio, leading to increased emissions of combustion products created due 
to the local oxygen deficit. A further consequence of the appearance of increased amounts of 
carbon monoxide (CO) and unburned hydrocarbons (HC) is to lower the combustion 
temperature, which determines the reduced NOx emissions. 

The deciding factor in the emissions of toxic compounds derived from transient states is 
primarily the value of force, which causes these conditions . But this is not the only factor. 
Another factor affecting the emission of toxic compounds derived from transients that has to 
be taken into consideration, is the condition of the engine. This condition, described with the 
structure parameters while using the engine, is constantly changing, which is responsible for 
the processes of wear. This change enhances the formation of toxic compounds during 
transient states, as these processes, though short, are so dynamic that the instantaneous 
concentrations frequently exceed ZT values of the steady states . Therefore, it is expected that 
the engine with its structure parameters changed due to wear, will be more sensitive to the 
effects of transients and thus it will be easier to determine its technical condition [5]. 

The basic parameter deciding about the correctness of combustion process in spark- 
ignition engines is the fuel injection timing. Even small deviations result in the significant 
changes in the key indicators of the engine, including the exhaust emission indicators. In the 
case of conventional engine design, a "self-acting" change of fuel injection timing is rather 
unlikely. However, in modern constructions, where most of the control parameters is 
controlled electronically, it 1s possible for the control system to be damaged and the settings 
of injection timing to be changed. 

The paper will present the modeling of transient states with a variable angle of injection 
timing and their impact on the changes in the basic concentration of toxic compounds. 


2. Identification of a dynamic process of multi-equation model 


Building on the experience of authors [6,7,8,9,10] with modeling of toxic compounds 
concentrations, it was decided to implement the multi-equation models, proven during steady 
state, for the analysis of dynamic processes, whereby it is assumed that the change process of 
gas toxicity occurs throughout a time, which makes it dynamic. Therefore, the model can be 
described as multi-equation system of linear differential equations. Since the measurement of 
the concentration of toxic compounds is a discrete measurement, discrete-time signal (time 
series) is a function whose domain is the church of integers. Thus, a discrete-time signal is a 
sequence of numbers. Such sequences are referred to as recorded in the functional notation. 
The adoption of such a notation was striving to minimize the impact of errors including the 
approximation of functions that would have to occur when using the continuous functions. 

Discrete-time signal x[k] is often determined by sampling x(t), a continuous signal in 


time. If the sampling is uniform, then x[k] = x(kT) . Constant T is called the sampling period. 


Course of the dynamic process in time depends not only on the value of force at a given time 
but also the value of extortion in the past. Thus, the dynamic process (system) has a memory 
where it stores consequences of past interactions. 
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The relations between the input signals x,[k],x,[k],---,x,[k], and output signals 
Yilk] yk]. У] , k=0,1,2,... , will be described by a system of linear differential 


equations. 


У[К *1]7 ay VLA] + аа, y, [X] + Вых + Вох [+ box, [X] + д, 


Inn 


ук+ 1] = а ук] + ау] ++ arn ук) + By X[k] bx; [k] +--+ + By, х, [k] + 6 (1) 


y,l[k *1] 2 a, VLD + a, y [+ a, y, [Kk] 5b, ox [A] + b, xs [K] + bpn lk] S, 


n mn n 


where: 


y,[k],i=1,2,...,m - output signal values at К, 

х [k], j 2 12,...,n - input signal values at k, 

a, — is a coefficient found in i-th equation with j-th output signal, i, j =1,2,---,m 

Б, - is a coefficient found in i-th equation with j-th input signal, i =1,2,---,m, 7 =0,1,---,7, 
Е, - is a non-observable random component in i-th equation. 


In analogy to (1), the system of equations (2) can be written in matrix form 


ylk + 1] = Ay[k] + Bx[k] +6 Q) 
where: 
b, by Dim a, dj Qn 
в b, by by, | А ау Any dy, | 
ba D e ба An а 
УА] plk +1] х [К] 9 
ves 9| sucio HE aol ы 
Vnlk] Jk +1] x,[K] Sm 
Later denoting: 
C := [A|B] = [cij]mxcm+n) (3) 


and 
= dul 
4d = [мар 
the system of equations (1) is shown in reduced form 


ylk + 1] = Cz[k] + € (4) 
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Identification of the system of equations (1) and (4) will be based on the selection of the 
coefficients using the set of measurements on the real object of input and output signals. The 
problem of aforementioned selection the authors present, among others, in [6,7,8,9,10]. 


3. Study of dynamic process in engine fuel supply system through multi-equation 
models 


The object of this research was the engine fuel supply system (fuel injection) of a single- 
cylinder test engine 1-SB installed in the Laboratory of the Exploitation of Marine Power 
Plants at the Naval Academy [8]. The experimental material was collected by trivalent 
developed a complete plan [4]. The implementation of specific measuring systems (measuring 
points) of the above experiment design were performed using a programmable controller, 
which allowed a high repeatability of dynamic processes. The period between an onset of the 
clipping of injection system components and the re-stabilization of output quantities was 
adopted as the duration of the dynamic process. This period was chosen through a series of 
experiments, and it averaged to about 106 seconds. 

In order to identify the impact of the technical condition of the fuel supply system on the 
parameters of the engine power during dynamic processes, sets of input quantities (preset 
parameters) and output quantities (observed parameters) were defined. For the purpose of this 
study a set of input quantities X was limited to three elements, that is: x; - engine speed n 
[r/min]; x» — engine torque 7, [N-m]; хз — fuel injection timing aww [ КЕС]. The study was 
conducted in accordance with the approved complete plan, for three values of speed, ie 850, 
950 and 1100 [r / min]. For each speed, torque (Ты) increased and thus created a transient 
state, consequently for the load of 10, 20, 30, 50, 70 [N]. For speed of 850 r / min, afraid of a 
large engine overload, the loads of 50 and 70 N were omitted. Similarly, this was done to the 
speed of 950 г / min and a 70 N load. Fuel Injection timing varied by + 5 °КЕС, yielding three 
values, i.e. face value - М, accelerated angle - W, delayed angle - P. 36 repetitive transients 
were obtained this way. Graphic interpretation of the test program is shown in Figure 1. 


—— 850 obr/min —— 950 obr/min ——11000br/min 


30 


t4 t2 t3 t4 t5 te t [s] 


Figure 1. Schematic of the program of research 


Similar treatment was applied to the set Y of output quantities, limiting the number of its 
elements to only the primary toxic compounds in exhaust manifold: y; — concentration of 
carbon monoxide in the exhaust manifold Ccow [ppm]; y2— concentration of hydrocarbons in 
the exhaust manifold Chow [ppm]; ys — concentration of nitrogen oxides in the exhaust 
manifold Суох [ppm], ул — tsp exhaust gas temperature [° C], ys— air-fuel ratio A. 
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Statistical identification was made using GRETL [1]. Estimation of the equation 
coefficients for specific output variables was performed using the least-squares method and it 
had to verify the significance of its parameters and, consequently, the rejection of 
insignificant values, which consequently led to a significant simplification of the models. 
Given the large amount of experimental material, for the purpose of this study 15 selected the 
most typical cases, while limited to the greatest loads that occur during the experiment. Tables 
1, 2, 3, 4 are coefficients of equations for nominal output variable injection timing. 
Significance values of model parameters are shown in the last column of the tables. Equations 
describing the changes in concentration of hydrocarbons (y2) and the concentration of nitrogen oxides 
(y3) have undergone the greatest simplification. (Table 1, 2, 3). In the case of equations describing 
changes in carbon monoxide and hydrocarbons, they depend on the excess air ratio A (у5) 
significantly, a parameter directly related to the parameter of the structure, which was the fuel 
injection timing (x3). Both CO and HC significantly dependent on each other. Furthermore, in the case 
of CO, speed has a greater impact, while in the case of HC it is the load, which seems to be logical, 
considering the creation processes of these compounds in the cylinder. 


Table 1. Least-squares estimation of the dependent variable y, 


Variable Coefficient Mean error Student t p value i significance 
xl I 0,785212 0,297116 2,6428 0,00952 pe 
y5 1 -640,111 84,9526 -7,5349 <0,00001  |*** 
y2 1 6,00076 0,536841 11,1779 «0,00001 |*** 


Table 2. Least-squares estimation of the dependent variable y; 


Variable Coefficient Mean error Student t p value i significance 
y5 1 52,7105 9,26814 5,6873 «0,00001 |*** 
x2 1 2,23258 0,472489 4,7252 «0,00001 |*** 
yl I 0,0605765 0,0102887 5,8876 «0,00001  |*** 


Table 3. Least-squares estimation of the dependent variable уз 


Variable Coefficient Mean error Student t p value i significance 
x2 1 0,104014 0,0251446 4,1366 0,00007  |*** 
xl I 0,216752 0,00166072 130,5167 <0,00001 |*** 
y4 1 0,025725 0,00437145 5,8848 «0,00001 |*** 


Table 4. Least-squares estimation of the dependent variable y; 


Variable Coefficient Mean error Student t p value i significance 
xl I -0,483745 0,23243 -2,0813 0,03994 |** 
yl I -0,0155244 0,00653628 -2,3751 0,01943  ** 


109 


y3 1 4,6947 0,975664 4,8118 «0,00001 |*** 


y5 1 -118,888 10,9896 -10,8182 «0,00001  |*** 


The graphical presentation (Fig. 2) of the matching can prove having a good model fit to 
the values obtained in the experiment on the engine, as well as equal distribution of residuals 
from the regression of mean values (Fig. 3). 
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Figure 2. Graphical matching of the CO model to the empirical data for the transient at n = 1100 r / min and 
load change with Ttq = 50 to N = 70 N Ttq a) delayed (P) injection timing angle, b) accelerated (W) injection 
timing angle 


The analysis of Figure 3 shows that the worse adjustment, due to the higher residue 
values, exists for accelerated (earlier) fuel injection timing angle. 

The results of the presented analysis highlight the significant advantage of multi-equation 
models, the possibility of multi-criteria analysis of the variables in the case where these values 
are in mutual correlation. Analysis of these relationships in one model reflects the reality 
more accurately (because there are obvious interactions between, for example, CO and HC 
and, for example, A), and thus allows for a broader interpretation of the test problem. 

Despite the undeniable advantages of multi-equation models do not provide direct 
information on the quality of changes, in this case, differences between the concentration of 
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various toxic compounds due to changes in fuel injection timing. Only the juxtaposition of 
courses of the experiment or the analysis of the obtained models provides some picture of the 
phenomenon. Nonetheless, analysis still remains difficult due to the similarity of transient 
waveforms, irrespective of the value of extortion. 
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Figure 3. Graph of the CO regression residuals for the transient at n — 1100 r / min and load change from 
Та = 50 N to Ти = 70 Ма) delayed (P) injection timing, b) accelerated (W) injection timing. 


In such case, it is desired to apply the criteria that would be useful in the objective 
assessment of the comparative levels or emissions from transients. The use of an evaluation 
index is one of the commonly used methods in such cases. The basic hourly evaluation is an 
individual emission of toxic fumes, which is calculated by the formula [3]: 


Е,=а, С, Gy: [gh] (5) 


where: 

j - СО, HC, NO,, 

а, — characteristic factor for a given compound j: 
асо = 0,000966, анс=0,000478, амох=0,001587, 

C,, — concentration of individual compounds [ppm], 
G,, , — exhaust gas flow [kg/h]. 


spi 
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This evaluation is, however, difficult to apply in the case of the analysis of transient 
processes, since determining the exhaust flow would require the estimation, thus introducing 
significant errors, which obviously excludes this method. Another view was proposed by the 
authors in [3], using the following relationship evaluation: 


t 

W, =a, |C (dt (6) 
0 

where: 

C;;(t) — the concentration of any toxic compound in time t [ppm], 

t — the duration of the transient state [s]. 


Thus, by integrating the area under the curve obtained from an experiment or from a 
model, an indicator that accurately describes the direction of change was obtained. On the 
other hand, this indicator continues to not describe the nature of the changes. As is known 
from observation, depending on the value of force, the course of the transient can vary 
significantly. These differences depend largely on the intensity of the experience of individual 
phases of the transient. Frequently, in the course of a typical transient, two phases can be 
noticed. First, characterized by the highest growth rate, accompanied by a sharp increase in 
the concentration of TS, typically several times more than the concentration in the steady 
state. The second phase of the transient is characterized by a much less violent course, it has a 
monotonic character and approaches the steady-state of concentrations in an asymptotic 
manner. 

As mentioned above, the concentrations of individual toxic compounds derived from 
transients are characterized by a certain regularity and repetition, and therefore a tool had to 
be found that would be deprived of the above-mentioned disadvantages of the indicators, 
while being able to be described in the precise and objective nature of the changes in the 
concentrations of individual toxic. It seems that the described method would be the analysis of 
the correlation of individual transients. This method determines the correlation of the 
researched transient state and that of the transient adopted as a model describing the 
phenomenon. Analysis of the correlation function allows you to specify the degree of 
correlation and its nature. Analyzing the components of the function can infer the said 
transient nature, that is, the participation and intensity of the individual phases. 

The graphic imagining of the analysis correlation is a scatter diagram presented in Figure 
4 and 5. Figure 4 shows the linear correlation function of the concentration of unburnt 
hydrocarbons HC at early injection timing angle ( 30 ? HVAC - red ) relative to the nominal 
injection timing (26 ° HVAC) , where the correlation coefficient was г = 0.75. The color green 
indicates the correlation function and the HC concentration at the delayed angle (22 ° HVAC) 
also with respect to the nominal injection timing. The correlation coefficient in this case was 
smaller, and was r — 0.59. Smaller values of the correlation coefficient were affected by the 
dispersion of points around the correlation function, which indicates an unstable transient 
process (matching multi-equation model is significant even in this case, as the greatest value 
of the residue is 60 ppm). Analogously, correlation analysis may be performed for NOx (Fig. 
5). 
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Figure 4. The concentration of hydrocarbons HC for the transient at n = 1100 r/min and load change from Ти 
= 30 Мю Ти = 50 М: P - late, W - accelerated injection timing, CHC (N, S, P) - HC concentration for (N) 
nominal, (W) accelerated, (P) delayed injection timing 
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Figure 5. The concentration of NOx for the transientat n = 1100 r / min and load change from Ти = 30 Мю Т, = 
50 N: P - late, in - accelerated injection timing, CNOx (N, S, P) - NOx concentration (N) nominal, (W) 
accelerated, (P) delayed injection timing 
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The correlation coefficients indicate a high matching of correlation functions in the two 
cases. And, for the accelerated angle (red), correlation coefficient is r = 0.79, while for 
delayed angle, г = 0.80. Another regularity can be noted, namely, much higher concentrations 
of NOx fall for the former angle, which of course is consistent with the theory of combustion 
(with the increase of injection timing, the share of kinetic combustion increases, at the same 
time increasing the pressure and temperature of combustion) , are thus improved conditions , 
at which nitrogen oxides are formed. For HC the situation is reversed (Fig. 4), delayed 
injection timing results in a higher proportion of diffusion combustion, which often goes into 
burnout and thus creates favorable conditions for the formation of unburned hydrocarbons. 


4. Summary: 


In the course of this study, the following conclusions have raised: 


e  multi-equation models provide a good match to the empirical results, 

e using a multi-equation model makes it possible to predict, and thus greaten the modeling 
of concentration change (emission) during the transient, 

e а discussion of the accuracy of different methods to estimate emissions, even using the 
method of spline functions is possible. 
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